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Abstract 


Several shell deformation models are developed for use in crashworthiness analysis of rotationally 
symmetric structures. These models use analytical techniques to predict the crushing force 
versus axial crush distance characteristics of both a rigid-plastic. hemispherical shell and an 
elastic. cylindrical shell loaded axially by a rigid flat plate. Additional methods are proposed to 
determine the effects of cutout sections and internal stiffening members on the crushing force 
capacity of the shells. These proposed methods are applied to determine the energy absorption 
capability of the composite/metal nose structure of a mini submarine subjected to a head-on 
impact with a flat rigid wall. The nose structure is composed of a stiffened composite shell that 
covers and is attached to the metal. forward, hemispherical portion of the pressure hull. The 
complex structure is simplified to a rotationally symmetric shell model and modified to account 
for the effects of stiffening elements and cutout sections. Laminated plate theory. a progressive 
composite failure method. and the models developed are then used to determine the structure s 
energy absorption capacity. The fairing structure is determined to be capable of absorbing the 
kinetic energy associated with an impact from an initial vehicle speed of 7.5 knots. Further. the 
entire nose structure (the fairing and forward pressure hull) is predicted to be able to absorb 
the kinetic energy from a 26.4 knot vehicle impact. 
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Сһартег 1 


Introduction 


The Advanced Seal Delivery System (ASDS) (Figure (1-1)) is a highly sophisticated. mini- 
submarine 19.8 m long that displaces 60 metric tons submerged. It is designed to transport a 
team of Navy SEALs from a host submarine to a tactical point of insertion and back. While 
the transported SEALs are a formidable offense weapon, the ASDS itself has no offensive or 
defensive capabilities and must rely on stealth and its ability to "hide" in the water column 
to ensure the safe transport of its personnel. As such, the ASDS will be operating in close 
proximity to the ocean bottom in order to mask its underwater signature. This near bottom 
operating strategy presents some problems. however. The most pressing of these is the increased 
risk of impact with submerged objects or structures. For the safety of the embarked crew, it is 
paramount that a determination be made as to the crashworthiness of the ASDS and whether 
operational restrictions need to be imposed. In effect. a safety envelope needs to be developed 
which takes into account the depth, bottom type and impact angle to give a maximum safe 
operating speed. A natural starting point for this analysis is to look at how a front end impact 
will affect the forward structure of the vessel. This analysis. however, presents no small task. 
The complexity of this problem quickly becomes apparent when the structure types and problem 


variables are viewed: 


e The structure of the nose section consists of various composite materials and steel in-. 
cluding: glass woven cloth. carbon cloth, kevlar. vinylester, rubber. syntactic foam, and 


HY-80 steel. 
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Figure 1-1: External Isometric View of the ASDS 
e The variables of the problem include: depth. impact angle, speed, and bottom type. 


Additionally, the interaction of the structural material and modes of failure during defor- 
mation and fracture are not well understood or predictable to a high degree of accuracy. In 
light of these complexities, using model testing to develop some empirical means of relating 
the structural behavior during a collision to key parameters such as impact angle, velocity, or 
crush distance would be the optimal analysis tool. However, due to the cost associated with 
such tests. this method is not feasible. Moving to the other end of the spectrum, a computer 
simulated "crash" test of the ASDS could, in theory. be used to ascertain the desired safety 
envelope at a substantially cheaper price and in an accelerated time frame. While modeling 
the full ASDS via computer simulation seems plausible on the surface, a little digging reveals 
that even the most sophisticated computer tools have not evolved far enough to accurately 
predict impact responses for such complex, composite/metal structures. With this in mind, it 
is apparent that neither model testing nor computer simulation alone can solve this problem. 
A middle road between these methods may be the required path. In effect, a combination of 
structural simplification and computer and analytical modeling will be used to complete the 
analysis. With all of this said, the overall objective of the project will be to determine a safe 


impact speed for the vehicle such that during an impact the pressure hull will not rupture and 





the crew will be able to survive the collision. 

As a first step in solving the problem, this thesis will analyze the energy absorbing potential 
of the composite nose fairing and forward pressure hull hemi-head structures during a head-on 
impact with a rigid wall structure. While this impact scenario is idealized and likely never 
to occur during actual operations. it will serve as a worst case impact and best case energy 
absorption scenario. This will be a worst case impact from the stand point that all kinetic energy 
will be absorbed by the structure in the form of deformation. No energy will be absorbed by 
the impacted structure or be maintained as vehicle kinetic energy. This is a best case energy 
absorption scenario in that the impact will allow the full utilization of the bow structure to 
absorb energy before the pressure hull is impacted. The reasoning behind looking at such a 
case first. is to ascertain the maximum energy absorption potential of the composite structure. 


This analysis will progress in three phases: 


e Deformation model development. 
e Structural modeling of the ASDS nose structure. 
e Energy absorption analysis. 


The models developed include deformation relationships for a rigid-plastic, hemispherical 
and an elastic. cylindrical shell loaded axially by a rigid flat plate. A method is also included to 
approximate the reduced load carrying capacity caused by cutouts in the shell structure. Next, 
the composite shell of the nose fairing will be modeled using simplified, rotationally symmetric 
shells (i.e.- cvlinders, frustum, and hemispheres). Internal stiffening members are incorporated 
into the model by “smearing” the effective area of each stiffener across portions of the inner 
surface of the shell to increase the effective thickness of the structure. Cutouts in the shell for 
components, such as thruster doors, video equipment, and lights are also modeled. Finally. the 
deformation relationships developed will be applied to the corresponding sections of the model 
to produce a force verse displacement relationship of the structure as it is crushed axially by a 
rigid flat plate. The energy absorption capability of the entire structure will then be determined 


and translated into the vehicle speed needed to produce such an impact. 
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Chapter 2 


Development and Discussion of 


Analysis Tools 


2.1 Crushing Response of a Hemispherical Shell Between Rigid 
Plates 


2.1.1 Prior Research 


Much attention has been devoted to the study of the energy absorbing characteristics of thin- 
walled. metal, dome structures [1]-[7). The problem of large deformations of rigid-plastic, spher- 
ical shells being compressed between rigid plates was first studied by Updike [1]. А simplified 
yield condition and geometry was used to obtain a sequence of limit loads, thus producing a 
force-displacement relationship valid from a displacement of a few shell thicknesses to about 
1/10th the shell radius. Wierzbicki [4] used similar geometric and deformation assumptions 
to derive a force-deflection relationship based on the optimization of bending and membrane 
energies valid theoretically until the plate deflection reaches half the radius of the shell. Com- 
parison of [1] and [4] show that the force-displacement results predicted by each differ only 
slightly. Experimental testing of hemispherical shells was reported by Kitching [5] and Kin- 
head [5]. Kitching et al studied the quasi-static loading of hemispherical shells and observed 
experimentally three different stages of deformation for shells with radius to thickness ratios 


ranging from 36 to 460. Later, Kinkead performed experimental testing on shells with radius 
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Figure 2-1: Geometry of Hemispherical Deformation Model. 


to thickness values ranging from 8 to 32 and identified two vice three deformation stages. The 
loading characteristic observed by Kinhead were compared to values predicted by the methods 
of [1] and [4] and large differences were found. Gupta [7] developed another prediction method 
using a rigid plastic analysis which relies on two calibration constants that must be determined 
from curve fits of experimental data.. The method proposed by Gupta showed good correlation 
with experimental results. This analvsis develops a means of predicting the load-compression 
relationship of hemispherical shells using a rigid plastic analysis technique with an emphasis 
on creating a simple. analytical model whose accuracy is improved over that proposed in [1] 
and [4]. The model developed will utilize the same deformation model used by [1] and [4], but 
calculate the membrane energy using a simple strain averaging method. Additionally, the re- 
sulting relationship will be capable of analyzing material with differing compressive and tensile 


stresses. 


2.1.2 Deformation Process 
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The spherical shell of radius R will be loaded from the top by a moving rigid plate with 
the base of the shell pinned to a fixed rigid plate. As the shell is loaded axially, a spherical 
dimple forms which is an inversion of the original shell curvature as shown in Figure (2-1). In 
this deformed state, surface AB remains rigid and undeformed. At point B, a circumferential 
plastic hinge increases the curvature of the shell, followed by rigid body translation between B 
and D and a reversal of curvature by the plastic hinge at D. The shell is subsequently unloaded 
and rigid body translation results on surface DE. As the shell deforms, the hinge circles at B 
and D move progressively outward. Deformation continues in this manner until the inverted 


portion of the shell at point E contacts the fixed rigid plate. 


2.1.3 Deformation Model 


Geometry 


This analysis begins with a look at the various geometrical relationships that can be attained 
from the deformation model in Figure (2-1). First, the radius to the center of the toroidal 
surface (a) and the displacement of the moving rigid plate (6) can be expressed in terms of the 


shell radius (R) and deformation angle (0) as 


а = sind. and (2.1a) 


6 = R(l-cosd) (21р) 


By differentiating Equations (2.la) and (2.1b), the radial and axial velocities of the toroidal 


section and rigid plate become 


da do 


pe em —— 22 
3i R(cos 0) at and (22) 
ав ао 

ne А а 2. 
2 R (sin o) m (92 


By approximating cos g(t) using the first two terms of its Taylor expansion in Equation (2.1b) 
and rearranging. a simple equation relating the deformation angle to the rigid plate displacement 


can be deduced. 


25 
u. 24 
Ф R es 
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Turning now to the toroidal surface BCD. the average plastic hinge length (L) and toroidal 


surface area (4) can be calculated by the Equations (2.5a) and (2.5b). 


= 27a, and (2.52) 


А = Алтао (2.5Ъ) 


Additionally. the radial and rotational velocities of the plastic hinge at point B can be taken 


as 
do 
= — 9 
V ny (2.6) 
ад С 
do _ YV_Rao em 
dt T TRGE 


Finally. in order to simplify subsequent calculation, surface BCD will be approximated as a 


parabolic surface of the form 


B M 
2 твіп о 





т? + т(1 – соѕ д) (2.8) 
with the x and y axes oriented as shown in Figure (2.2.3). 


Derivation of Membrane Stresses and Fully Plastic Bending Moment 


Next. the membrane force (~) and fully plastic bending moment (M) relationships for 
a material having differing compressive (Ce) and tensile (o) flow stresses will be developed. 
As shown in Figure (2-2), when 0, Æ o+ the neutral axis will shift away from the geometrical 
center of the cross section by an amount 77. By using the fact that the neutral axis of a specimen 
subjected to a bending moment is located where the membrane force (.V) is equal to zero, we 


are able to calculate a general expression for 7 in terms of shell thickness (f). c. and o; as 


5+1 2-1 
NT ih o,dz — | сва» =) (2.9) 
0 JO 


where both c; and c, are taken to be positive.. Upon integrating and solving for 7, Equation 


shown below 
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Figure 2-2: Membrane Stress and Fully Plastic Bending Moment Representation for Material 
having Differing Compressive and Tensile Flow Stresses. 


(2.9) becomes 


1 6--С; 
= — 2.10 
7 2 0.+0+ | 


With 7 calculated, general expressions for N and M can be derived as follows: 


t 
№. = / ordz 
0 
2+1 n 
o,zdz + сс242 
0 0 


Note that the regions of membrane deformation are compressive in nature while the bending 


М 


forces are both tensile and compressive. After integrating, substituting Equation (2.10) for 7, 


and setting a — 2 D the total membrane force and fully plastic bending moments become 


N = Ot ДІЛ) 
l > ot? 
M = -0.4 —— = —a (2 12) 


and (2.12). o. can be set equal to o; which 


mt 


In order verify the validity of Equations (2.11 
results in a symmetric stress profile throughout the specimen. In such a case Equations (2.11) 


and (2.12) reduce to the standard expressions for the fully plastic membrane force and bending 


moment. 


= 
сл 








Energy Balance 


The principle of virtual velocities will be used to relate the rate at which energy is being 
introduced to the system to the rate at which it is being absorbed. Equation (2.13) states that 
the rate at which energy is applied to the shell must exactly match the rate at which it is being 


absorbed by the bending and membrane effects in the shell. 


1) (2) (3) 
dé dE, E ШЕТ 


— = oe o 
dia di wat (2-13) 


The remainder of this section will further develop each of the terms in Equation (2.13). Starting 
first with the left hand side of Equation (2.13). term (1) can be rewritten in a different form 


by using Equation (2.3) as follows 


dé do 
— = Iba аа 9 
PT = PRsino” (2.14) 


The second term of Equation (2.13) represents the rate at which energy is absorbed by the shell 
due to bending processes. Using the deformation model shown in Figure (2-1), it is plainly seen 
that bending energy will be absorbed at two discrete locations - plastic hinges B and D. This 


component can be quantified by the following expression 


dE» 40 
= p Mi La (2.15a) 


Equation (2.15a) can be further simplified by making the following assertions 


e the magnitudes of the bending moment and rate of angular rotation at hinges B and D 


are equal. 


e the total hinge length $` JL; can be rewritten as 2L (where L is the average hinge 
i=B,D 
length. 
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Figure 2-3: Deformation of Surface BCD. 


Using these simplifications and Equations (2.1a), (2.5a). and (2.7), Equation (2.15a) can be 


expressed in the following form: 


dE, Я с у М до 
ке үү ш. 9. 
T 47R* sing ms (2.16) 


The final term in Equation (2.13) accounts for the membrane energy dissipated on surface BCD 


and can be written as: 


ШЕ), ағ 
—— = NA— 27 
dt dt = 
d=: Е : ; SEDEM : . Е 
where © is the instantaneous rate of strain experienced by surface BCD. As a simplification, 
the average strain rate, Fan will be used in place of the instantaneous strain rate. Figure 
—— 
(2-3) shows a detailed view of the deforming surface BCD. The average strain of this surface 


will be calculated by dividing the surface into two subsurfaces, calculating the average strain 


on each. and then averaging the strains from each region to get the final average strain. This 








statement can be expressed mathematically as 


oy Slat SI 
Eau =, = re 


(2.18) 


— 


where ¢;,,, and €7;,, are the average strains experienced as surfaces BC, (region 1) and C,D, 
(region 2) deform to BCand CD respectively. Looking first at region I, the strain can be written 


in terms of initial and fina] curvatures of the BC. 


ах 


~ 


1 / ду u 
ТЕГЕНЕ) er (2.192) 


In general. the average strain in region I can be formulated as: 





1 rsin 
Ela = | стат (2.20) 
0 


rsinó, 


Upon substitution of Equation (2.19a), integration and simplification, Equation (2.20) reduces 


to 


Elav = 56? (2.21) 


The strain in region II requires an additional step to calculate. While the strain in region I 
is purely compressive for the entire deformation, the strain in region II is initially compressive 
for the first part of deformation and then tensile for the remainder of the deformation. So, 
the average strain will be calculated for the compressive phase as surface CD. is deformed to 
an intermediate flat surface CD and also for the tensile phase as the shell further deforms to 
surface CD. With this. the average strain in region 11 can be calculated in the same manner 


as in region 1. The details of the calculation follow. 








Ellas EE С ШЫН еј = Ша, ы 

ii rsinó 
Elav E T. sin о І (еен Е А) F [ere ш || = 

К 2 

1 rsin® I, 1/4 

EN = = So = 
rsino Jo 2 2 \ ах 

_ Е 2 2 (ә) Фу) 
| (2.22) 
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With the average strain calculated in both regions, we can now determine the average strain 
— 

of the entire surface BCD by substituting Equations (2.21) and (2.22) into Equation (2.18). 

Subsequently, the average strain rate can be calculated by simple differentiation of the average 


strain as follows. 


Гу 
Eav => p (20199) 
dE av do 
di 9d к 


Now a final expression for the rate at which energy is absorbed by membrane forces can be | 
realized by substituting Equations (2.5b), (2.1a) and (2.24) into Equation (2.17). 


ШЕ е же a 
——— =4лт\ e > 
m 4т"Х Есіп оо m (2.25) 


New expressions have now been formulated for each term in Equation (2.13). By substituting 
the new expressions (Equations (2.14), (2.16), and (2.25)) into Equation (2.13) and solving for 


the instantaneous crushing force (P), the following result is obtained 


2772 


r 


Рат 


Further. the mean crushing force (P) can be calculated and minimized with respect to the 


toroidal radius (r) to obtain an optimal toroidal radius. 








== T R М 2 
= = 2.97 
A s (5 +9) р 
аР NR 
"uA 2 (570 -1)-9 
1 N 
2 2.28 
Е- ?V RM 52 


By substituting the optimal toroidal radius back into Equation (2.27) and using the fact that 
at the optimal radius the membrane and bending contributions are equal, the mean crushing 


force can be expressed more simply as 
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e N 
5- = 99 
Popt = 40R\/ 5 (2.29) 


Equation (2.29) can be rewritten in terms of plate displacement and shell thickness by using 


Equations (2.4), (2.11), and (2.12). 
26 
= и 9) | 
орї 8 [л (2.30) 


йе | 


Finally. the optimal instantaneous crushing force becomes 


Ром = AnV/20 02 t262 223) 


2.2 Derivation of Cylindrical Tube Equation 


2.2.1 Prior Research 


The deformation and energy absorption characteristics of thin-walled. composite, cylindrical 
shells have been studies for many years by many researchers. [{8]-[15]. Farley {10]-[13] and Hull 
[14] have performed extensive testing on composite shell samples composed of various mate- 
rials using different lay-up conditions. Both quasi-static and dynamic loading situations were 
investigated. These studies did much to further the knowledge base in the areas of determining 
crushing modes and their controlling mechanisms. They also showed that there was no explicit 
set of rules that could be generalized for all composites types. Material lay-up pattern. crushing 
speed, fiber volume fraction, and radius to laminate thickness ratio all play key roles in deter- 
mining the deformation characteristics of the shell. From these studies. complicated numerical 
models have been established for use in predicting the failure of very simple composite struc- 
tures. The degree of success realized by these models has varied in even simple situations. The 
vast majority of the work described has been directed toward understanding the deformation 
of composites on the micro- and macroscopic scales. Hanefi and Wierzbicki [15] took a more 
global approach to a similar deformation problem. They developed a simple model that was 
used to predict the energy absorption characteristics of a combined metallic/composite tube 
structure in axial compression. The model coupled a rigid plastic analysis of the metal tube 


with a compressive perfectly plastic and tensile elastic analysis of the composite material. The 
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Figure 2-4: Deformation Model for Cylindrical Thin-walled Structure. 


predicted results matched well with experimental data. In this analysis, a global failure ap- 
proach will also be taken. A deformation model will be developed using elastic principles to 
determine the force deflection and energy deflection characteristics of a thin-walled, cylindrical 
shell loaded axially. Bending and in-plane stiffnesses will be step-wise varied to account for 
stiffness variations as plys fail within the structure. Finally, a simple result will be established 


for a material exhibiting constant and identical flexural and inplane stiffness properties. 


2.2.2. Deformation Process 


Figure 2-4 shows a representation of two buckling lobes of the cylinder during deformation 


where 6 denotes the axial crush distance, a is the radius of the cylinder, and 2H is the half 
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length of the buckling lobe generated. As the cylinder is crushed axially by a distance 6, the 
cylinder will form an axisymmetric, elastic buckling pattern with a wave length equal to 2H. 
Deformation takes place both in the form of bending and circumferential strain. The material 
will behave elastically until its elastic limits are reached with subsequent deformation occurring 
with the material in a progressively degraded state. The material degradation takes place 
through a step-wise reduction of the material’s flexural and circumferential elastic moduli and 


material thickness. The process continues progressively through all buckling lobes formed. 


9.2.3 Deformation Model 


The principle of virtual work (Equation (2.32)) will be used as the starting relationship from 
which to derive an expression relating the crushing force (P) to the axial crush distance (6). 
Equation (2.32) states that the energy applied to the cylinder through an external force (1) 
must be exactly equal to the energy absorbed by the structure in the form of circumferential 


(2) and bending (3) strains. 


(1) (2) (3) 
De || Т (2.32) 
р 
ов Вер апас = (2:953) 
Еа |. (егде -- Еггағ:даУ (2.34) 
. 


Because the deformation will be elastic, Equation (2.32) can be simplified using the elastic 
relationships between stress and strain (Equation (2.33)) to obtain Equation (2.34). Figure (2- 
5) shows a detailed view of the deformation of a single buckling lobe with parameters pertinent 
to the analysis labeled. Using this figure, expressions can be derived in terms of the variables 
shown for the incremental volume (dV), strains (zz and eg). and incremental strains (de; and 
deg). First. inextensibility will be assumed for the shell as it buckles, making dV constant over 


the deformation process and having the following value: 


dV = 4natH (2.35) 
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Figure 2-5: Detailed View of the Deformed Cylinder Section. 


23 








Next, the circumferential strain can be expressed as the ratio of the outward deflection of the 
shell (w) to the original shell radius. The actual deflection is measured from the original surface 
of the shell to the deformed surface and varies from 0 at the being of the lobe to a maximum 
value of w, at the midpoint, and then back to 0 at the end of the lobe. In order to simplify 
future calculations, a linear deflection profile (w(x ) = — is used to approximate the actual 
deflection profile as shown in Figure (2-5). Using this linear relationship. £g and deg can be 
expressed as 


, 


и T 
— — = ——- M ә 5 
29 5 que (2.36) 
т 
gg = ——4@4и„ 287 
а 0 zA и ( 31) 


The bending strain component can be determined by using the curvature of the shell section 
(x) and the position within the thickness of the material (z) as is shown in Equation (2.38). 
This variation of the bending strain throughout the shell thickness will be assumed linear with 
maximums occurring at 5 апа -$ and no strain occurring at the mid-plane. Further, plane 


sections will remain plane as shell curvature increases. 


(2.38) 


Er = ZRK = 


z 
R 
However, the radius of curvature (R) is related to the half lobe length and the maximum 


circumferential deflection by Equation (2.39) 


Ко (Ко је = H -u 
eo 8 (2.39) 





Substitution of Equation (2.39) into Equation (2.38) yields the final expressions for ¢, and 42, 


(Equations (2.40) апа (2.41)). 
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Er = ШО (2.40) 
2 
de, = ішо (2.41) 








Now. Equations (2.35)-(2.37) and (2.40) can be substituted back into Equation (2.34) to obtain 


the following: 





(4zat) E 3 
pat (ба ШЕ | Em(x )?dx + алан) f zaz) wodwo (2.42) 
0 


If the flexural and membrane Young’s moduli can be assumed equal and constant throughout 


the deformation process, Equation (2.42) can be further simplified to Equation (2.43 ). 


Рао Е 


ІН at? 
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2 wodwo (2.43) 


Previous work [18] has proven that the elastic, buckling lobe length of a cylinder can be related 


to the shell radius and thickness by 


2H = 1.72Vat 


0 
© 


Н (2.44) 


Using Equation (2.44), Equation (??) can be reduced to 








4. tat 
d == ) d ) 
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A final simplification can be made to Equation (2.45) by relating ô to wọ. Using Figure (2-5) 


the following geometric relationships can be generated. 


біп(о) = H (2.46) 
Соза = 1-27. (2.47) 








By squaring and adding Equations (2.46) and (2.47), the relationships are combined in the 


following form 
bo 6 


2 2 
1= =] +11-— 2 
q] + | 2 Еј (2.48) 
After differentiating, rearranging and substituting Equation (2.44), the desired relationship is 


achieved 


моди = va - 2 = 


This expression can now be substituted back into Equation (2.45) to obtain the final expression 


relating the crushing force applied to the axial deflection of the cylinder. 


4 
P(6) = ¿ME = WE 4 (2.49) 


As previously states, Equation (2.49) applies to materials with equal and constant flexural and 
circumferential Young's moduli. If these conditions are not met, the more general expression 


(Equation (2.42)) must be used. 


2.2.4 Correction for Open Sections 


Open sections or cutouts in a shell will effectively reduce the load required to deform the shell 
in the local region of the cutout. These cutouts could even change the local deformation mode 
of the structure or change the order in which the overall structure deforms. For simplicity. this 
analysis will include only the reduction in loading capacity for several reasons. The relationships 
previously developed. Equations (2.26) and (2.49), are valid only for the failure modes assumed: 
so. a change in failure modes (i.e.- from axisymmetric shell yielding to lobar buckling for ex- 
ample) will not be captured here in. While this seems like a fairly unreasonable assumption. 
the reader is reminded that the objective of this analysis is to determine the mazimum energy 
absorption capability of the structure. The development of Equations (2.31) and (2.49) assume 
failure modes which produce this maximum. So. a shift to a different mode would produce a 
lower energy absorption. With the assumption that the failure mode remains unchanged at the . 
cutouts. it does not make sense to consider the order in which parts of the structure collapse 


because the same amount of energy will be absorbed regardless of the collapse sequence. So. the 


26 








only effect that will be considered is the reduction of strength caused by the loss of material. 
This effect will be approximated by multiplying the applied load of an intact shell by a reduc- 
tion factor in the local vicinity of the cutout. This reduction factor will be proportional to the 
remaining shell area divided by the material present in an intact shell of the same dimensions 


(Equation (2.50)) 


ТО лола = БІЗ нестане (0) where 


(6) = К а 


2.00 
A complete shell (б) 5 ( ) 


While this area ratio is part of the reduction effect, it does not capture all aspects. Other 
factors. such as cutout geometry. affect the load carrying capacity of the structure around 
the cutout. While these other effects are recognizable, a theoretical derivation of their effect is 
impractical. This effect must be determined experimentally. The experimental correction factor 
would be used to match experimental test data with theoretical predictions and is expressed 
as € in Equation (2.50). Since no experimental calibration was done during the course of this 


thesis. € will be set equal to 1 until future evaluation determines its correct value. 


Hemispherical Dome 


For the hemispherical dome. the cutout sections will be assumed to be circular with diameter 
D as shown in Figure (2-6). From Figure (2-6) it is apparent that the correction factor can be 


stated in terms of the average plastic hinge radius (a(6)) and angle 0(6) as follows: 


age. erago Hine. Dimni r Enei ee Onen en 
p Average Hinge Diameter hemisphere 
2ла(6) - Ф(б)а(6) o(6) 
=(6) zu Bm hemisphere =|1- Әт. Ehemssphere (2.51) 


The task now becomes to express o(6) in terms of known quantities. Figures (2-7) and (2-8) 
show more detailed views of the side and top of the shell respectively with pertinent dimensions 


labeled in each figure. 
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Figure 2-6: Top View of Hemispherical Dome with a Single Circular Cutout Section. 
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Figure 2-7: Geometrical Relationship of Parameters Used in the Derivation of the Correction 
Factor for Open Sections in a Dome Structure (side view). 


From these figures. the following geometric relationships can be established: 


Al = rg-- Dcos — a(6) 





A1 
2 
= сов0 
АЗ = Р-А? 
Р D 
2 
Ase (5) _ да? (2.52) 


o(6) can now be expressed in terms of the initial radius at which the cutout is encountered (го). 
cutout diameter (D), cutout inclination angle (@), and average plastic hinge radius as shown in 


Equation (2.53). 
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Figure 2-8: Geometrical Relationship of Parameters Used in the Derivation of the Correction 
Factor for Open Sections in a Dome Structure (top view). 
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(0) - 2arcsin (22) 


| ( то+ совв—а(6) fet Dees sate) | (р _ (EE) 
сов 0 cos 8 (2.53) 


oló) = 2arcsin a6) 


Next. a(6) can be related directly to the radius of the hemisphere (R) and the crush distance 


of the shell (6) by the following: 


alô) = /R? - (R б)? = VSR- ô) (2.54) 


Finally. substitution of Equation (2.54) into (2.53) and the result into (2.51) yields the desired 


result for the correction factor (Equation (2.55)) 


то + D cos0—4/6(2R—6) D то + 0 созд6—,4/6(21—6) 
cos O = cos O 


=(6) 1 : arcsin Е 
— — — 1 ин mas r 
т . v&(2R - 6) wi 


(2.55) 


Cylindrical Shell 


The same principles described in Section (2.2.4) will also be used to estimate the effect of 
cutout sections in the cylindrical section with the exception that the cutout shape will be much 
simpler. The cutouts in the cylindrical section will be rectangular in shape, and thus, the 
material removed will not vary with axial crush distance. So, using Equation (2.50) and again 


assuming that € = 1. the correction factor of the cylinder can be expressed as 


27a — noa n 

----- -1- Е where (2.56) 
n = Number of cutouts 
> = Angular extent of cutout 
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Figure 2-9: Stiffener Smearing Process. 


2.2.5 Effect of Shell Stiffeners 


The energy dissipation effects of shell stiffeners can be taken into account using Equations (2.31) 
and (2.49) by assuming that these structural members are “smeared” across the surface of the 
shell so as to produce a new effective shell thickness. This process is illustrated in Figure (2-9). 
If the material properties of the stiffeners are identical to those of the shell material. their cross 
sectional area is distributed across the surface of the shell such that the new effective thickness 


(terr) can be determined by the following expression 


n 
> А 
tarr = == 2.57 
ff = tshet + = (2.57) 


where A, is the cross sectional area of the ¿t? stiffener. n is the total number of stiffeners. a is 
the shell radius. and tsre is the shell thickness. If the shell and stiffener materials have different 
properties. the force distribution over the effective cross section must be made consistent with 


the force distribution prior to the “smearing” process by using the following relationship 


Ре} = Pheil + Pstif feners (2.58) 
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By substituting the appropriate relationships and assuming that Pers will be calculated using 


the shell properties, Equation (2.58) can be written as 


п 
27at,r | Езћепе = гпаћ ће Езћене + > АЁ е 


тс 


which can be solved to determine the effective shell thickness 


to 
сл 
e 
— 


1 n 
t z — A E, А 
eff = tshet + 5— Е 2- : ( 


When E; — E» —...— E, in Equation (2.59), Equation (2.57) is recovered. 
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Chapter 3 


Modeling and Analysis of ASDS 
Composite Shell Structure 


With the analysis tools developed. focus now shifts to the modeling and analysis phase of 
the project. First. a model will be developed for the actual ASDS nose fairing structure and 
forward hemi-head of the pressure hull using simple, rotationally symmetric, geometric shapes 
for the shell. Stiffeners and cutout sections will be incorporated using the techniques discussed 
in Section (2). Next. the analysis tools developed in Section (2) will be applied to determine 
the crushing force to axial crush distance relationship and energy absorption capability of the 
structure. Finally. a maximum impact speed will be calculated by equating the energy absorbed 


in the collision to the kinetic energy possessed by the vehicle just prior to the impact. 


3.1 ASDS Modeling 


3.1.1 Actual Structure 


The forward portion of the ASDS consists of a hemispherical pressure hull and a composite 
fairing structure. The hemispherical pressure hull section is constructed of HY-80 steel and 
is welded to the cylindrical portion of the hull to form a water tight boundary. The ASDS 
nose fairing (see Figure (3-1)) is a composite structure attached to the forward end of the 


pressure hull by eight titanium clevises. The structure has an overall length of 2.94 m and a 
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Figure 3-1: The Actual ASDS Nose Fairing Structure Showing Internal Components (Viewed 
from Aft Looking Forward). 
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Table 3.1: Properties of Materials Used in ASDS Nose Fairng. 


pseudo-rectangular cross section with maximum dimensions of 2.21 m high and 2.06 m wide. 
The structure itself consists of an outer shell constructed from E-Glass (1583) /Epoxy (7780) 
preimpregnated (prepreg) woven cloth in a single mold process. Internally, the shell is stiffened 
by a network of longitudinal and transverse woven carbon cloth (6K 5HS Satin Weave IM?7) 
beams which are bolted to the shell and interconnected by titanium clevis fittings. Attached to 
the front end of the shell is a sonar dome comprised of a BFGoodrich Rho-Cor material system. 
This system consists of Kevlar (K285)/Epoxy (7714) prepreg woven cloth with a rubber core 
material. The properties for the various material systems described are contained in Table 
(3.1). Cutouts have been made in shell and the sonar dome to accommodate such equipment 


as thrusters, the forward anchor, lights and video cameras. 


3.1.2 Simplified Model 


The actual structural design of the ASDS had to be simplified to rotationally symmetric 
shapes before the analysis tools developed in Section (2) could be applied. A representation 
of the model developed is shown in Figure (3-2). For the forward pressure hull, no modeling 


idealizations needed to be made since the structure was already rotationally symmetric and 
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Figure 3-2: Model of ASDS Forward Section (Dimensions in meters). 


hemispherical in shape. The pressure hull has a .99 m radius and an undisclosed thickness. The 
actual thickness of the structure is held as confidential by the U.S. Navy so as not to reveal 
the operational diving depth of the vessel. For this analysis. a notional thickness of 19 mm 
will be used. The sonar dome was modeled as a spherical dome section using the axial depth 
of the actual sonar dome and a base radius equal to the average of the height and width of 
the actual sonar dome. The thickness and lay-up construction was taken to be that of the 
actual structure. Three cutout sections simulating the locations of the camera and lights were 
located on the model 120? apart at the same radial distance. The actual and model sonar 
dome dimensions are summarized in Table (3.2). The shell structure was modeled using a 
combination of a frustum for the forward portion of the shell that transitions into the sonar 
dome and a stiffened cylinder for the after section. The radius of the cylindrical section was 
taken to be the average of the maximum height and width of the actual structure. The frustum 
was used as a transition section between the sonar dome and the cylindrical section and uses 
the corresponding radius at each end to fix its dimensions. For the frustum section. the shell 
thickness was taken to be that of the actual material. The thickness of the cylindrical section 
was taken to be an effective thickness calculated using Equation (2.59) to account for the added 


structural stiffness of the longitudinal carbon beams. The transverse beams were not utilized 
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Table 3.2: Comparison of Actual and Modeled Sonar Dome Dimensions. 


in the effective thickness calculation because their energy absorption in an orthogonal impact 
is would be minimal. Instead, these beams were used to define nodal points for buckling waves. 
Table (3.3) summarizes the actual and modeled dimensions of the shell section. Two cutout 
sections 1.13 m high by 
openings created by the forward thruster doors. (It is assumed that during a collision imminent 
situation. the thruster doors would be open and the thrusters in use.) Additionally, the titanium 


clevises (used to attach the internal beams together and to the pressure hull) and the anchor 
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.81 m wide were included in the cylindrical section to account for the 


were treated as rigid objects. 


3.2 Analysis 


3.2.1 Impact sequence 


As stated previously in Section (1). this analysis will investigate an head-on impact with a. 
flat. rigid wall. Structural deformation will initiate at the bow as the sonar dome is impacted 


first. The deformation will progress .635 m through the sonar dome, .95 m through the frustum 


38 













A ы | 
[306 | Width ај БЕНЕН БЕНЕН 
A Average [ — 
Г Г ваш 1107 
35 


Frustral en 
| Гепа (т) | ћ ( 


Fwd Radius ( 

ee [9 
Aft Radius ( 

[| соды веде авы | ^^ 


Table 3.3: Comparison of Actual and Modeled Nose Shell Section Parameters. 























section of the nose shell, and .18 m into the cylindrical section before the anchor makes contact 
with the wall. At this point, some of the crushing load is transferred to the pressure hull via 
the anchor and simultaneous deformation of the cylindrical section and the pressure hull results 
for another .66 m. The deformation process ends when the after most transverse beams and 


titanium cleaves make contact with the wall. 


3.2.2 Sonar Dome 


Equation (2.31) will be used to determine the crushing force required to deform the structure 
throughout the progressive crushing of the sonar dome. Implicit in the use of this equation is 
that the structure can be modeled using rigid-plastic techniques. While it is understood that in 
general composite structures behave in a more brittle than ductile way. kevlar composites have 
been shown to buckle under axial loading in modes similar to those seen in metal structures 
119). This behavior coupled with the fact that the structure has a rubber core was used to 
justify the use of this analysis technique on this component. 

The sonar dome is composed of 22 plys of kevlar cloth laid-up symmetrically sandwiching a 
3.81 mm rubber core ((90. 0];. 0, rubber core, 0, [90, 0]5) for a total laminate thickness of 9.4 mm. 
The compressive and tensile flow stresses were calculated using a combination of laminated 


plate theory and a progressive failure analysis [21]. The detailed calculations are shown in 
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Appendix (A) and, for brevity sake. only an summary of the process and the result will be 
included here. First, the material properties of the ply (listed in Table (3.1)) were used to 
determine the undamaged 0° ply stiffness characteristics. These properties were then modified 
using techniques of [21] to determine the reduced stiffness of the 0° ply with matrix and fiber 


damage. The stiffnesses were calculated using the following relationship: 


E; и Е} 
l—viv2 1—viv2 


Ze E» E 
Q- [3925 тт (3.1) 


0 0 Ез 


The stiffnesses were then rotated to a 90? plv orientation. The stiffness of the rubber core 
material was also determined from Equation (3.1) by setting E] = Ez. vı = va. and E3 = IL 


The overall laminate stiffness was then calculated using the following: 


| a. | 
А = — — (Nos hos Qoe + Море ћорг 90° + Мецфегћтибђег Стићђет) Ућеге 


Пашниме 
Nz = Total number of “x” plys 
= Thickness of individual ply 


Q = Stiffness of individual ply 


Next. stress space strength parameters were calculated from the material test data provided in 


Table (3.1) using Equations (3.2) and (3.3) and the generalized von Mises value of —3 for F D 


1 Fra 
хх V X X'YY! 
uA Fla EM 2 
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КЕН ЛА (3.3) 
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The strain space strength parameters and failure envelopes could then be determined through 
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the relationships in Equations (3.4)-(3.6). 


Сы = FiiQuiQi (3.4) 
(E T Ge = 1 (3.6) 


A progressive failure model [22] was then used to approximate the stress strain relationship of 
the material under pure compressive and tensile loadings (oriented in the 0° ply direction). The 
progressive failure model allows for a rough prediction of post first ply failure behavior of the 
laminate. The model works by calculating the strain of first ply failure for the assumed load 
condition and then determines the stress using the undamaged laminate stiffness ({a} = A*{<}). 
After first ply failure, the limiting ply is failed to a degraded state, a new laminate stiffness is 
calculated. and the next ply failure strain and stress is determined. A ply can be damaged at 
most twice — once by degraded matrix and the final by fiber failure. In certain cases. based on 
the properties of the plys used, the degraded ply failure strains can be less than that of the 
undamaged ply. For this condition. the ply failure occurs after its undamaged failure envelope 
is surpassed. The tensile and compressive stress versus. strain estimates are shown in Figure 
(3-3). 

The compressive (7,5) and tensile (Cto) flow stresses were then determined to be 8.41 MPa 


and 53.35 MPa respectively using Equation (3.7). 


1 es 
Ef Jo 


Now. the crushing force (Psonar dame) and cumulative energy absorption (Esonar dome) OÉ the 
sonar dome can be determined from Equation (2.31) and (3.8) with the cutout correction factor 


(Equation (2.55)) applied to the crushing force over the appropriate regions of the shell. 


é 
E= | P(6)dé (3.8) 


The resulting crushing force and cumulative energy absorption profiles for an axial crush dis- 


tance (6) are shown in Figure (3-4). Here it is seen that the total energy absorbed for a crush 
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Figure 3-3: Compressive and Tensile Stress vs Strain Estimates Calculated using the Progressive 
Failure Method. 


distance of .635 m is .058 MJ. This figure also clearly shows the effect of the cutouts on the 


crushing force curve. 


3.2.3 Frustum Section of Nose Shell 


The frustum section is composed of a 22 ply. E-Glass (1583)/Epoxy (7780), woven cloth com- 
posite system laid up symmetrically ([90, 45, 0, —45]5, 90, 45, 0, 0, 45, 90, [90, 45, 0, —45]2) for a 
laminate thickness of 10.1 mm. The analysis of the frustum section was done by approximating 
the frustum as a cylinder with a radius equal to the average of the forward and after radii of the 
frustum and using the analysis discussed in Section (2.2.3). The details of these calculations are 
contained in Appendix (B). In order to use this technique, the appropriate circumferential (Eg) 
and flexural (Ef) Young's moduli needed to be determined. The method described in Section 
(3.2.2) was initially used to determine the strain space failure envelopes for each of the ply ori- 
entations used in the laminate. Subsequently. the failure envelopes were used to develop an E» 
versus. tensile. transverse strain relationship using the progressive failure technique described 


previously. E»? was chosen as the circumferential strain Young's moduli because as the buckling 
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Figure 3-4: Crushing Force and Cumulative Energy Absorption as a Function of Axial Dis- 
placement for the Sonar Dome. 


lobes are formed and bulge outward and inward, the laminate will be strained in the 2-direction 
with respect to the laminate orientation. In reality, both compressive (bulging inward) and 
tensile (bulging outward) strains will be experienced and thus, different E> versus. strain re- 
lationships will result. However. because the compressive and tensile strengths of the plys are 
approximately equal (Table (3.1)), the difference is small and only the tensile E» versus. strain 
relationships will be used to simplify the analysis. Figure (3-5) shows the E» versus. strain 
relationship that was determined as well as the order of ply failure. 

The actual circumferential strain was then calculated as a function of x and 6 using Equa- 
tions (2.36), (2.48). and (2.44). This relationship is shown in Equation (3.9). 


т тб Vat 


1 
с 3.9 
РЕ a тес 52 
Since the circumferential strain varies both with x and 6 and E2 varies with 6. Eg will vary 
with both z and 6. Eg was then calculated for discrete values of r and 6 in preparation for a 


numerical integration scheme to be used in Equation (2.42). 


To determine Ey, a relationship between flexural strain (ef) and the axial crushing distance 
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Figure 3-5: Transverse Young’s Modulus Variations with Pure Tensile. Transverse Strain for 
Frustum Section. 


(6) was first developed using Equations (2.38), (2.39), and (2.48). (see Equation (3.10)) 


2 
e= Eva- (3.10) 


Equation (3.10) was then used to determine a ef versus. 6 curve for each ply over the entire 
crushing distance for one buckling lobe (2\/at). For each. the strain was determined at the 
midplane of the individual plys. Using these plots and the strain failure envelopes, the dis- 
placement at which each ply failed in tension could be determined. Figure (3-6) shows how this 
process was done for the 0° oriented plys in the laminate. 

Flexural stiffness (D*), Ey. and laminate thickness (t) were then recalculated after each ply 


failure using Equations (3.11)-(3.13). 


D [figa (3.11) 


Ер = а (3.12) 
t = nh where (Sum 
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Figure 3-6: Flexural Strain to Axial Crush Distance Relationships for the 0 Degree Oriented 
Plys in the Frustum Laminate. 


п = number of intact plys 


h - ply thickness 


Figure (3-7) shows the variations in Ey and ¢ as 6 increases. Notice also that as more plys fail, 
E; becomes more erratic and diverges further and further from the in-plane modulus. This 
is expected since the laminate starts out homogenized and gradually de-homogenizes as each 
successive ply fails. 

The crushing load and cumulative energy absorption can now be determined as a function 
of axial crushing distance (Figure (3-8)) using Equations (2.42) and (3.8). Since neither Eg 
nor Ey was constant or equal during the deformation process, the simplified form of Equation 
(2.42) could not be used here. Figure (3-8) shows that five buckling lobes form as the frustum 


is crushed and that a total energy absorption of .105 MJ is realized. 
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Figure 3-7: Flexural Young’s Modulus and Laminate Thickness Variation with Axial Crush 
Distance of the Frustum Section. 
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Figure 3-8: Crushing Load and Cumulative Energy Absorption Characteristics of the Frustum 


Section. 
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3.2.4 Cylindrical Section of Nose Shell 


This outer shell of the stiffened. cylindrical section is composed a 28 ply. E-Glass (1583)/Epoxy 
(7780), woven cloth composite system laid up symmetrically 
((90. 45, 0, —45]2. 90, 45,0, 0, 45, 90, [90, 45, 0, —45]2) for a laminate thickness of 12.8 mm. Ad- 
ditional support is provided by an internal. frame structure constructed using carbon. woven 
cloth (6K 5HS satin weave IMT), composite, C-channel beams. The beams are 32 plys thick 
for a laminate thickness of 11 mm and laid-up symmetrically in the following configuration: 
[90.45.0.-4514[90.45.0.-45]4. Тһе shell analysis was accomplished by the same process as de- 
scribed in Section (3.2.3) and is shown in detail in Appendix (C). However. the effect of the 
stiffeners still needed to be taken into account. This was done by “smearing” the stiffeners 
into the shell and then calculating an effective thickness (tess) using Equation (2.57). The 
detailed stiffener and tes calculations are included in Appendices (C) and (D) with the teff 
determined to be 35 mm. An effective buckling lobe length was then determined based on teff. 
The crushing force and cumulative energy absorption could then be determined using Equations 
(2.42) and (3.8). Figure (3-9) shows the force and energy profiles at various displacements and 
also. that the cylinder was only allowed to crush .84 m in two buckling lobes while Table (3.3) 
indicates that the section is 1.36 m long and should buckle in three lobes. This reduction in 
length was due to the effects of the transverse beam structures. The forward beam and titanium 
fittings served to force a nodal point at the beginning of the thruster door section as seen in 
Figure (3-2). The after beams and rigid titanium clevises at the aft portion of the thruster door 
create another node point and serve to terminate the deformation process. So, only the .84 m 
shell And stiffener portion between the forward and after beams will deform during the impact. 


.49 М.) of energy are absorbed during this deformation. 


3.2.5 Analysis of Forward Hemi-head of Pressure Hull 


The forward pressure hull hemi-head is a metallic structure constructed from HY-80. As pre- 
viously stated. the thickness used in the analysis will be assumed to be 19 mm due to the 
confidentiality of the actual thickness. The deformation starts as the forward anchor contacts 
the collision boundary and is pushed back into the pressure hull. The deformation is assumed 


to be of the form described in Section (2.1) and, as such. Equation (2.31) will be used to de- 
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Figure 3-9: Crushing Force and Cumulative Energy Absorption of the Cylinder Section of the 
Nose Shell. 


termine the crushing force. The flow stress was calculated from a partial stress versus. strain 
curve using Equation (3.7). The detailed calculations and the stress versus. strain profile used 
are included in Appendix (E), and the crushing force and cumulative energy absorption versus. 
displacement results are shown in Figure (3-10). Throughout this analysis. it is assumed that 
the forward hemi-head remains water tight. Further analysis must be performed to determine 
the point at which the pressure hull fails. By allowing the pressure hull to deform to the full 


axial extent allowed for the nose shell, 7.566 MJ of energy will be absorbed. 


3.2.6 Kinetic Energy Calculations 


After determining the energy absorption capability of each modeled nose fairing component, a 
determination of the kinetic energy possessed by the ASDS was needed for various speeds. The 
kinetic energy was determined using Equation (3.14) 


P 1 
KEasps = x(Mvehicte + added) U^ (3.14) 


һә 


where Muehrcle 15 the submerged mass, Mgddeq is the added mass, and v is the speed of the 


vehicle. 
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Figure 3-10: Crushing Force and Cumulative Energy Absorption for the Forward Pressure Hull 
Hemi-head. 
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Figure 3-11: Description of Standard Submarine Weight Terms 
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Parameter Weight 
Pome | mete tons) 


Table 3.4: ASDS Weights 






The submerged mass of the vehicle was calculated using the ASDS weight report [16]. 
Weight reports for submarines are standardized and broken down into categories of weights 
based on the operating condition of the submarine. Figure (3-11) shows a brief summary of 
the standard terminology used in submarine weight analysis and how each relates to the other. 
Condition A is the dry vehicle weight. When the variable load (variable ballast tank water. 
personnel, etc.) is added, the new weight created is the normal surfaced condition (NSC). The 
main ballast tanks are filled with water to submerge the vehicle: so, the weight of this water 
added to the NSC becomes the submerged displacement (4,45). Finally, during submergence. 
non-water tight free flood areas (i.e.-between the fairings and the pressure hull) fill with water. 
This water. while in pressure equilibrium with the surrounding water. moves with the submarine 
as it transits, and when its weight is added to Asub the resultant is the envelope displacement 
(Aeny). For this case, it is this Aen, that is required for the analysis. Table (3.4) shows the 
various weights used to determine the envelope displacement. 

Next. the added mass was determined. Since this vehicle is relatively slender, it is expected 
that the added mass in surge will not be large. However, two approximations were made using 
a sphere and an ellipsoid. The size of the sphere was set such that the frontal area was equal 
to that of the ASDS (Figure (3-12)). Once the sphere dimensions are known, the added mass 


can be calculated using Equation (3.15) 


1 2 


с ЧЕ 3 Ee e 
Madded sphere = 5 Pwater Vsphere = 3 TPwater Rsphere = 4.5 metric tons (3.15) 


metric tons 
(where Pwater = 1.025 z) 


The added mass of the vehicle was also approximated using an ellipsoid with the frontal area 
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Figure 3-12: Dimensions Used During Added Mass Calculations 


equal to that of the ASDS (as in the case of the sphere) and a major axis length equal to that 
of the overall vehicle (Figure (3-12)). The added mass of the ellipsoid can then be determined 


using Equation (3.16) and the graphs from [17] for the added mass coefficient (Ca) 
4 a | | 
ed ellipsoid ~ Caz TPwaterab = 2.1 metric tons (where Cr, ——10) (3.16) 


Since the actual vehicle shape is not as bluff as a sphere. yet not as fine as an ellipse, the added 
mass of the vehicle lies somewhere between the two calculated values. So. the average value 
was used as the vehicle added mass. Now, the kinetic energy possessed by the ASDS at various 


speeds can be determined using Equation (3.14) as shown in Figure (3-13). 


3.3 Results 


With the energy absorption of each section determined. the overall energy absorption of the nose 


structure can be determined by summing the energy absorbed by each section throughout the 
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Figure 3-13: Vehicle Kinetic Energy for Various Speeds. 


impact sequence. Figures (3-14) and (3-15) show the integrated crushing force and cumulative 
energy absorption profiles for the entire deformation process. Two curves are shown in Figure 
(3-15) depicting the energy absorption of the nose structure with and without the effects of the 
pressure hull included. Additionally. the kinetic energy curve shown in Figure (3-15) is overlaid 
on this plot. Using the curves in Figure (3-15) initial speed versus. crush distance assessments 
can be made. Without the pressure hull included. the nose structure can absorb .67 MJ of 
energv which would allow the complete absorption of all kinetic energy possessed by the vehicle 
with an initial speed of 7.5 knots. If the effects of the pressure hull are included, a total of 
8.235 MJ of energy can be absorbed. This equate to an initial vehicle velocity of 26.4 knots. 
As previously stated, the point at which the pressure hull will fail during its deformation was 
not determined: so. a the safe impact speed for the vehicle cannot vet be determined. However. 


when this point is determined. Figure (3-15) will provide this safe speed estimate. 
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Chapter 4 


Summary and Conclusions 


Throughout the course of this analysis. several new deformation models have been developed. 
The force required to deform a hemispherical shell loaded axially with a rigid plate by axisym- 
metric dishing was related to shell thickness, material flow stress (compressive and tensile), 
and axial deflection using a rigid-plastic analysis technique. Additionally, a deformation model 
for a cylindrical shell composed of elastic material. crushed axially by a rigid flat plate, was 
proposed. This model assumed a concertina type deformation mode and used shell thickness, 
cvlinder radius, material stiffness, and axial displacement to determine the force required to 
crush the shell. Further, methods were set forth to include the effects of cutout sections and 
internal stiffening members on the crushing force characteristics of the shells. The effect of the 
cutout section was taken to be a reduction of the intact shell crushing force by a ratio equal to 
the shell area remaining divided by the total area of the intact shell in the local vicinity of the 
cutout. The effect of internal stiffeners was incorporated by increasing the thickness of the shell 
material to a new effective thickness determined by making the load carried by the effective 
structure equal to the combined load on the stiffeners and shell. 

These tools were used to determine the crushing force and energy absorption characteristics 
of the composite nose fairing structure and forward pressure hull of the U.S. Navy’s Advanced 
Seal Delivery System in a head-on impact with a rigid wall. The complex structure of the 
ASDS was first simplified and modeled to permit the use of the tools described above. The 
shell structures were approximated as rotationally symmetric shells (i.e.- cylinders, frustum, 


hemispherical domes) with the internal shell stiffeners being incorporated into the model in 
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two ways. The energy absorption of the longitudinal elements was captured by increasing the 
effective thickness of the shell in the vicinity of the longitudinal beams, while the transverse 
stiffening elements served only to force nodal points in the shell buckling pattern since their 
energy absorption during a head-on collision is minimal. Further, major internal components. 
such as the titanium clevises (used to connect the internal beams) and the forward anchor. were 
treated as rigid. massless objects. 

The structure was divided into a hemispherical sonar dome, a frustum and stiffened cylin- 
drical nose shell, and a hemispherical forward pressure hull based on the shape. material com- 
position. and thickness of the actual structure. Using laminated plate theory and a progressive 
failure model. stiffness versus crush distance and material thickness versus crush distance rela- 
tionships were developed for each of the material systems used. These relationships were then 
used in conjunction with the rigid-plastic hemispherical and elastic cylindrical shell models 
to predict the force deflection and energy absorption of each component. The component- 
level energy absorptions were then summed to approximate these characteristics for the overall 
structure. Finally. the absorbed energy was equated to the vehicle's pre-impact kinetic energy 
in order to determine an initial speed just prior to the collision. The composite nose fairing 
structure was deemed capable of absorbing the impact energy associated with a vehicle speed 
of 7.5 knots. If the effects of the pressure hull were included. the structure could absorb the 


kinetic energy from a 26.4 knot impact. 


4.1 Areas Additional for Study 


This analysis has taken the first step in determining the energy absorption capability of the 
ASDS. However, considerable work remains to be done before a final assessment can be made 


of the vehicle's overall crashworthiness. The following is a list of topics for further study in this 


project: 


e Model validation and calibration: The models developed herein have not been cali- 
brated or confirmed via actual test data. Cylindrical and hemispherical composite sections 
of the appropriate length to radius and thickness to radius ratios need to be constructed 


for each composite material and crushed under quasi-static and dynamic conditions. The 
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results should be used to confirm. improve. and calibrate the methods developed in this 


thesis. 


Rupture Point of Forward Pressure Hull: The overall objective of the project is to 
determine a safe impact speed for the vehicle such that during an impact the pressure hull 
will not rupture and the crew will survive the collision. The rupture or collapse point of 


the forward pressure hull needs to be determined and may vary with depth of the vehicle. 


Strength Analysis of Nose Fairing Mounting Plates: The composite nose fairing 
is mounted to reinforced metal plates that are welded to the forward hemi-head of the 
pressure hull. An analysis needs to be done to determine whether these structures will be 


capable of withstanding the force of the impact or if they will bend or shear off completely. 


Energy Absorption of Internal Components: The energy absorption potential of the 
internal components (such as clevises, angle bulkheads, the anchor, the sonar equipment, 
etc.) have not been included in this analysis. The effects of each of these components 
can be determined individually and added to the curves developed in this report to get a 


more accurate approximation of the energy absorption of the vehicle. 


Effect of Impact Angle on Energy Absorption: This analysis looked at an “impact 
sequence” for a head-on impact angle. In reality. this situation is likely never to happen. 
A relationship should be developed for various impact angles using model testing or other 
means to relate the reduction in energy absorption that is realized by an off-axis impact 


to that of an orthogonal impact. 
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Appendix A 


Sonar Dome Analysis 
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Ply Stiffness and Layup Inputs: 


1. Material System: Kevlar (K285)/ epoxy (7714) woven cloth 


2. Layup Configuration: [90,0]5,0,0,[90,0]<. 


Dig: 12 N 245-0 

М 90 :=10 N p45:70 
3. Total Number of Plys Total plys:7 N o4-N n454- N p45+N 90 
4. Ply Thickness һ:2.0і-іп h rubber; 15 in 
5. Laminate Thickness t [ 77 Total plvsh 4-h rubber 


іт 9.3989mm 


6. Strain Applied to the Material (circumferential): 


B=atan(e,/e,) 8 :=0 deg 
Shear Strain: E ¿:=0 


7. Number and Orientation of Ply Angles: 


Piys:=4 
0 
a:=1..Plys ‚_| 90 
ОЕ ЕЕ | ес 
:=1..3-Plys и 
c-=1..3-Ply _45 
8. Engineering Constants: 
Kevlar Rubber 
E ,:74 7.MSI Е; := 45.М$1 E rubber 71304 psi 
Eq:734GP :=.499 
MS шу ~ rubber 
G m :=1.26-GPa 
у оо m 
Ут.=5 vo 
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Model Parameters 


1. Radius of Base of Sonar Dome: 
2. Axial Depth: 


3. Radius of Curvature 


Radius to Center of Cutouts 


> 


5. Axial Distance to Center of Cutout 


о 


. Slope of Dome Surface at Cutout Center 


7. Diameter of Cutout 
8. Radius where Cutout Starts 
9. Radius where Cutout Ends 


10 Number of Cutouts 


Ply Strength Parameter Inputs: 


R model :=36. ап 
D model :=25-1n 


2 2 
._1 D modei +R model 


R curvature шы poem R curvature 7 0.976*m 


a” 
R cutout/ 727 in 


| 2 2 
D cutout *=R curvature— a/R curvature — R cutout D cutout = 0.282*m 


өст -D cutout R curvature 


A D cutout N D cutout t2- R curvature 


o =58.003%deg 


г 0 :=23.6934п 


гр. =29.664п 


.... 
D cutouts ' 7? 


1. Tensile Strength in Longitudinal Direction Х:=56.К$1 
2. Compressive Strength in Longitudinal Direction ТІ 
^ prime * ; 
3. Tensile Strength in Transverse Direction 
Ү :=56.К1 
4. Compressive Strength in the Transverse Direction 
Y :=14K 
5. Shear Strength рпте 714 RSI 
S :24 5-KSI 


Matrix and Fiber Damage Inputs 


БЕ 52-0 n y:=5161 n:=1 


E m_star = 15 n s:=316 


Е ху star: 5 
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Damage Ply Calculations: 


E s damaged F) 


v о 
ШЕП Еа 185: 
Ga -= = 
a Es Om 
р -] 
“т 
ny 
M SV 
Vm Er F Em 
= 2 
Vim = 
n P 
da ee 
“т ОКР Ст 
EE g= 
vf 


O [ma Ply Data: 
ndamaged ply, ply with matrix damage and ply with fiber damage. 


Ex 


Ex‘ 


Ex 


EXE [ star 


Vx 


V xx^ 7| V x E m star 


E. 


E vy ‘=| Ey damaged E m star 


E y, damaged E f star 


F xy. star 


FF xy. star 7| F xy, starF m star 


Vx E f star F xv. star f. star 
YY, ј=у ҮҮ prime ү ргіте SS, jS 
Rotation Matrices: 
ml, ‘cos 9, n, sin 9, 
e 2 
т! п. 2ml,n, 
P 2 2 Э 
Ir n4 ml, -2ml,n, 
2 2 
-ml.n, ml,n, т! = 
Rotation! а = 
^ 2 
ml, ~ i 
a 
Rotation2, ‘= n, = mi, 2 
mln, -ml n 


Es 
E ss :7| E s damaged E. m star 
E s damaged E f star 


XX 


SX 
q 


XX SX С 
prime, рите’ Е г 





e Evo 
M . Vives qe 
Ууд Ет 
4 
4 4 А 2 2 
n, 2- ml, Sn, 
9 
E ml E IN 2 
a 
2 2 2 “| 4 
п. т] n, min: 
2 l 2 2 = ] 2 2 
na m n4 == mi n. 
Pm -ml -n 3 ml -n 3_ т! 
а а 
3 3 3 
a - ml т, т! п -ml n 


Ут 
er 
Em 

2 2 

4 ml n, 

2 2 

4 т n = 

а 

2 2 

-4 ml, uS 








Ply Stiffness Calculation: 


Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 


Қ Е У -E 

E хха V yyq E xa у rubber = rubber une 0 

1-у ЛТ 1-У хы УУ | — rubber 1—У rubber. 
A E ._| Y rubber E rubber E rubber 

ООЛ “ху буу, Ууд ; Q rubber | ————- 0 

[= a m Ууд те Ey Ууд 1—¥ rubber 1— V rubber 
0 0 E 0 0 2E rubber 
554 IA 
1+Y rubber 


Creates the stiffness matrices for every combination of ply angle and ply state. Qnew is a 3 x 4 matrix 


9а ІШІ 
Өле Trotg,a, 9042; 9гца; 
Q rot, „:=Коапоп1._. ха = 
ara 5 M^ 12 © пема 2 Q года а 3 Q rot 2 О rotg à 6 
Q, T Q rotg a Е Q rotg 4 р O rota a П 


Stress Space Strength Parameters: 
separate F must be caiculated for each ply state 


1 FF Xy, Starg A 
XX XX 2 
д^ рите xx, XX prime, YY q YY prime, 
ER FF xy_Slarg 1 
: а | ee a ana eS YYYY e 2 
А IX, кх primeg VY, YY primeg qo 
1 
0 0 
E 
SS 
q 
1 1 
хх XX primeg 
Fy = 
4 1 


| 


Yy YE 
q primeg 
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Strain Space Strength Parameters: 


1. Strain strength parameter for each ply state 


8 
a iil % 1,1] | 
Ес 2 
2 9 81,2 ДЕ 
а? 1 1 аа | 
* Fas so] 99 1] “ға, 
G baselg ' 3a | D Sors 
2 Де 
COP em бат 90 2+[ А i 
* Fa 22 9% 9, 9% 2.2 а 
0 
Е Е 
G NEN Сатр ва 5 9 > 
Баѕе2 ' 
Шы четы баз > 


, 
, 


] Q 
a 


t3 


q 
ч 


il 


F M 
gue dq І 94 i 


2. Calculating G matrices for all combinations of ply state and angle: 


G Базе! | 
С basel l 
Rotation! 


базе! | 


С Базе! | 


G Базе 15 
G basel > 


1 | 
Rotation], : 
G basel 2 


G basel; 
G базе 15 


Rotation! 
basel3 


G basel} 


G Базе! | 


G basel» T 
25 


G basel 3 3 


to 


to 
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to 


G 2trans'” 


Rotation2 


Rotation2 


Rotation2 


1 pj 
Rotation2, 


Rotation2, 


а: 
Rotation2, 


1 Ў 
Rotation2, 


Rotation2, - 


кә 


to 


t 


кә 
ии 


кә 





G ltrans ‘> 


Rotation |, : 


Когапоп 1 


Rotation S 


Rotation! 3 


Rotationl, 
> 


Rotation 1; 


Rotation] 4 


Rotation |, 


Rotation |, 


1 


G basel] 


G basel] 3 


С basel> 
О разе1ђ 


G basel> 


G basel» 
= 3 


С pasel5 
G basel 3 
G basel 3 
С разе1 3 
G basel] 
G basel] 
G basel] 
G basel] 
G basel> 
С Ьаѕе15 


G basel> 


С basel> 3 


G basel3 
G basel3 
С basel 3 
G basel3 
С вазе! | 
G basel] 


G basel] 


С вазе! 1 3 


G basel> 
С basel> 
С pasel, 
G basel> 
G basel; 
С разе1 3 
С basel 3 


G basel 3 


— 


Когапоп2,. 
С base23 32 
с base2] 1 

Rotation2, : 
Ба5е21 , 
G base2> 1 

Rotation2, - 
Бае22 , 
Обоз | 

Rotation2, · 
base23 , 
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бал | Giwang, У ша. Уа) 
О = G Itrans, 4 C Itang 2 G lans. € G 2trans. 2 


б шап. ; Gltrany ¿ У ша) Уа, 


Formation of Strain Space Failure Envelopes: 


1. The following relationship is used to determine the failure envelope: 6; pret Gee! 


2. After shifting to polar coordinates, the above equation can be solved for r using the quadratic equation: 


> 
Е |r cos £5=F-sind Ат + Вт 120 


3. Calculation of coefficients for the quadratic formula 


A(6.6):= G, , зе)“ +2: 6, y 3050) sin(9)+ Gp, уох) 


B 6,b.g, :-2 б, _SIN(O)-E + б, А SEs G, i ЗЕ st G, | gem) 
2 ~, , , 


~ 
= 


3 
о": | 





_-B 6.b.e5 «AB 6.0.2, ~—4-A(6,b) C Ф.Б; 
Es TACo.b) 


Strain Space Failure Envelopes: 





6 :—0-deg,l.deg..360.deg v(x) ‘=tan( B)-x х.-0, 


Ply Strain Strength Envelopes 


Transverse Strain 








Longitudinal Strain 


——— 0 plv-Undamaged 
- O plv-matrix damaged 
-----. 90 plv-Undamaged 
---- 90 ply-matrix damaged 
: 45 ply-Undamaged 
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Calculation of Degraded Laminate Stiffness: 
e first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 


stiffness calculated based on the ply damage history. 


plot! := jale | 
ble! 
count 2 
nl] 
se 1 
while 1 
el ¡0 


ol jt 0 MPa 


£l oun min | 


ty B,al,es 
Еу В.Ы+3,Е; 











count! count ant B) 
l В rubber 

A Sou Ваље [x 0° Qneway 1 5+3 90° © пем er Q rubber | 

Mos No03 

=] -1 
E Мош | ^ státcount | т | 
El count 
Ol count EOS Dunt count 
al—al+2 if el oun €x 8.41. 5 
blebl+!1 if ne х В,Р1+3,Е 5 
counte count- l 
break if а1-- 0156 
break if al>3 
break if bl>3 
tempstoree- augment el 91 
8 Р "MPa 
Е ‚_[-236 © rubber; a 
ТЕҢБЕ? 1 li -=1.. rOWS € rubber © rubber;;; ‘=| Q rubber 0 plot:=stack plotl, augment € = 
0 І 
Tensile Flow Stress Determination «х. :=пһер роќ!>,роҝ2> ‚x 
| <1> 
| plot rows plot) 
ЕЕ. curve( x)d x.MPa 
| <і> 0 
plot rows plot) 


а | = 53.354°MPa 








Calculation of Degraded Laminate Stiffness ЕЕ 
e first ply failure strain 15 first calculated. Subsequently, a progressive failure of plys is allowed and the 


stiffness calculated based on the ply damage history. 


В :2180-4ес 
plot] -= |а1< 1 

ble] 

counte- 2 

nlel 

sel 

while | 
ol j € 0 MPa 
acount, 1-4! 
а оит, 2 P! 


1 mi 
E count 1^ 


E x Balies 
ex ß.bl+3,e, 








count El count an B) 
1 ћ rubber 
А Staf-gunt Бр E 0° Q пема] | $+М№ 90: 9 пем о SF ncs Q rubber | 
er , , 
NOS 0S q 
=i =l 
E паши || A statcount |, | 
£l count 
2 
бі оша” А Stafcount count 
aleal+2 if em Ex B,al,e, 


blebl +1 if ele еу В,Ы1--3,66 


counte count 4- 1 
break if al+bl=6 
break if al>3 
break if bl>3 


ol 
tempstoree- augment £l, 





MPa 
06 шп 9 rubber 
Е al | | iii .=1.. rows Е rubber б rubber; = О rubber: 0 plot :=stack plot!,augment є rubber —yqp3— 
0 | 


Compressive Flow Stress Determination 


<l> 
plot rows plot) 
. curve( x)d x-MPa 
0 
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x 
curva X) :—linterp plow! > plot ^? EX 


ЗЕ 1 
сот 


u с с= 8 419МРа 
plot rows plot) 





Force vs Axial Crush Distance Determination: 


1. Establishing Displacement Array .—D model { bbb— 1) 


bbb' nn— | 


2. Calculate Open Section Correction Factor 
AS) =N: 2R curvature — 9 


T o d cutour Cos e ) — al ô ) 
a CE ШОШ 


ao) 


T od-d cutout €0s( 6 ) — a(6 ) 


E an очо) 


N cutouts 9(0) 


=(0).=1 — >= 


3. Calculate Unbalance Flow Stress Correction Factor 


Er. 
T а = 1728 
Zero: 


4. Calculate Crushing Force for Displacement Array 


1 
2 


1 
Е. бы < otn Ao eg: бырь 


= 
~ 


кој [97] 
tu — 


if бы S pe ô bbb К са 
Pl(xx) :=(Iınterp &.P,xx)) 


Energy Absorption of the Structure: 


Епегем( 611) :=итмегр 6 ,Р,611) 


5 
Сшпшапуе Елегеу,,, = Energy( x)dx Cumulative_Ener, ows) = 0.058°MJ 
0-т 
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пп:=30 _ Б6р:=].. пп 
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"ENS. 
. O bbb 44/26 са 
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Appendix B 


Frustum Section Analysis 
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Ply Stiffness and Layup Inputs: 


1. Material System: E-Glass (1583)/Epoxy (7780) woven cloth 


2. Layup Configuration: [90,45,0,-45],,90,45,0,0,45,90,[90,45,0,-45], 


М0:56 N p45 i74 
N 90:56 N pas :=6 
3. Total Number of Plys Total_plys:= Мо-+М 145+ М pas+N oo 
4. Ply Thickness h:7 0184n 
5. Laminate Thickness t [ ;7Total plysh 
ty -00101%п 


6. Strain Applied to the Material (circumferential): 


P=atan(e,/e,) В :=89.5 ее 
Shear Strain: E ¿:=0 
7. Number and Orientation of Ply Angles: 


Plys :=4 


0 
a:=1.. Plys 90 
Б 6:=| _ |-deg 
45 
c:=1..3-Plys А 
-45 
8. Engineering Constants: 
E , '=3-MSI E ¿:=.55.MSI G m :=1.26:GPa 
E m :=3 4GPa 
E, :=2 7.MSI m 
У pem b 
“Vm 5 
uc 


ЇЇ 


sequence .= 


deg 








Component Parameters: 


1. Radius of Shell: 
2. Length of Shell: 


3. Wave Length of Elastic Buckling Wave: 


4. Number of Buckling Waves Along Shell 


Ply Strength Parameter Inputs: 


1. Tensile Strength in Longitudinal Direction 


. Compressive Strength in Longitudinal Direction 


. Tensile Strength in Transverse Direction 


2 
3 
4. Compressive Strength in the Transverse Direction 
5 


. Shear Strength 


R shell ,—3253.ft 


L shell ЗА" 


Н: = shet L 


.—L shell 
Number waves ‘= TH 





L shell 


Number waves *=> Не: 


X :=39 1-KSI 
X prime /733.8 KSI 
Y :=39.1-KSI 
У рите :=33.8-К51 


$ .=9 41.К $1 


Matrix and Fiber Damage Inputs 


Е f star :=01 Ny :=.5161 


Ет жаг" 15 145-5316 


Damage Ply Calculations: 


m 
^ 


Ут 
| m 
E v. damaged = EI Ix Е 
т ; m 
PERI | 
n "m 
> 
А | 1 1 ME 
E се Е).т === 
s damaged F) um GUT m 
ae 


ul 


n:=1 








F xv. star: 7-2 


Number waves = 4 733 


2-Number waves 











Ш Ply Data: 
ndamaged ply, ply with matrix damage and ply with fiber damage. 


| E. ESSE 


EXE f star 


v N 
ухх У ХЕ m_star 


Vy Ef star 


ҮҮ Y YY prime, Y prime 


Rotation Matrices: 


ml , -=cos En n, :7sin 9, 
= ә 
т! n, 2ml,n, 
з= 2 m? -2ml 
а п ma = Тата 
2 
-ті п. ті п ml — п 
а а а а 
т! = n 
a 
n 4 ml " 
a 
2 2 2 
ml n, ml -en 
Rotation], = 
ml 2. п 2 ml 2. п 
а а 
3 
т! т -ml -n 
ml .n - ml 
a a 
2 2) 
ml, * n, * 
a 
Rotation2, :— 2 2 
= п 
а а mi, 


ES 


E y. damaged E m star 


E y. damaged E f star 


F xy. star 


FF xy. star? 7| F xy. starE m. star 


F xy_starE f_star 


554 :25 
2 

2- ml п 
2 

2 ті ПА 


Es 
E ss := E s_damaged E m_star 


E s damaged E f star 


XX .=Х 


XX SX Е 559 
1 ш E 
pr теа prim E Р 





Ем 

v Ууд: хет 
Аза 

4 ml 2, ^ 

4 ml 2 NS 

2 2 
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Ply Stiffness Calculation: 


Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 


E хха у Wa Е x: 
|— у = V Yq l-v y Ууд 
9. :=| Уха У E] е 
l=v ха У js l—v = Уу Ууд 
0 0 Е Sq 


Creates the stiffness matrices for every combination of ply angle and ply state. Qnew is a 3 x 4 matrix 


Q i1 
Q 
Q, Кр Гоа | Q rot, a 3 ne 5 
Q rot, ., ;" Rotation] - B a 
q,a a E = О пека a Cor: 3 Q rotg a > Ca 6 
Q Q Q 
E 33 г014,а 5 [Olga 6 О, а 4 


Flexural Stiffness: 
1. Indexing variables ff:=1..Total_plys+ ! 
gg .=1.. Total_plvs 


== 


2. Establishing an array of thickness positions for each ply 





= + 


3. Assigns the appropriate stiffness matrix based on the layup sequence 


Qfex,,'7 O new; 1 if sequence P0 deg 


g 
if =90.d 
і sequence, ер 


55 
9 new) 5 
О new] 3 if sequence, 545 deg 


Q new| 4 otherwise 
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Stress Space Strength Parameters: 
separate F must be calculated for each ply state 


FF 
Xy, Stato 


А | XX XX primey YY q тт prime, 


= FF Xy_Starq | 
a У 0 
ХХ КА prime YY ҮҮ primeg q prımeg 


| 
XX XX prime, 





0 0 > 
SS, 


Strain Space Strength Parameters: 


1. Strain strength parameter for each ply state 


i Fag 25 412 


E 
ое 


С зал 
base2q pe 0 E Br E 


Sud hee q 4 22 
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2. Calculating G matrices for all combinations of ply state and angle: 


G Itrans ‘= 


Rotation . 


Rotation] 1: 


Rotation | ii 


Rotation l4 


Rotation l,- 


Rotation l,- 


= 


Rotation |, - 
> 


| 


G basel | 
С вазе! | 
G basel] 
С вазе! | 
G basel> 
G basel> 
G basel» 
G basel» 
G basel3 
G basel; 
G basel 
G basel3 
С Базе! | 
G basel] 
С вазе! | 
С вазе! | 
G basel> 
С һазе12 
G basel> 
G basel> 
G basel; 
G basel3 
G baselz 
G basel; 
G basel] 
G basel | 
С вазе! | 


G базе! 
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С 2trans = 


G Баѕе2 

Rotation2, - 
u G Баѕе2 
G base2> 

Rotation2, : 
G base2> 
G base23 

Rotation2, . 
G base23 
G base2 | 

Rotation2, - 
чк Баѕе2 
G base2> 

Rotation2, - 
21 G pase2, 
G base2; 

Rotation2, - 
21 G base23 
G base2 

Rotation2, - 
3 G base2] 
G base25 

Rotation2, 
G base2> 
G base23 

Rotation2, - 
3 G base2 3 
G base2| 

Rotation2 y 
G base2 | 
G base2> 

Rotation2, - 
G base2> 
G base23 

Rotation2 y 
G base23 








Rotation! 
G basel> 


G basel> 


с 
G basel 3 L] 
G basel3 E 


J 
G basel^ 11 
G basel> А 


Когапоп l5 
G basel 3 


G baselz А 
С вазе! | 
G basel | ~ 5 
Rotation], 22H 


G базе 1 | 


G basel l 


Rotation! a 
G basel^ 


G basel> 


G base; 11 
С ает; 555 


Rotation! 4 


G baselz 12 


G basel 3 Р 


G basel> L] 
G basel> E 


G шапа | G lrans ; G шап; ; G 2trans. | 
Gee G Itrans. 3 G Itrans. ә G Itrans. 6 G Ztrans. 2 


С lrans 5 С гап. 6 С Itrans. 4 G 2trans, 3 








Formation of Strain Space Failure Envelopes: 


1. The following relationship is used to determine the failure envelope: G; О! 


2. After shifting to polar coordinates, the above equation can be solved for r using the quadratic equation: 


3 
erco 5-5г51пф A^ 4B r— 120 


3. Calculation of coefficients for the quadratic formula 


2 . 2 
A(6,b):7 G, Ego a с, | о ши G, ү? 


Во БЕ. 2-2: G5 Е БӘЛЕ G5 Е um б, i O с, | wa 


de y 


5 
Cop. E. «= Soc EE с = |] 


1 


Bobe, +в bi AG ee, 


2-А(9.5) 


Strain Space Failure Envelopes: 


6 :70-deg,l.deg..360deg v(x):-tan(B)x Ce cog) £x 6,b,g£, ‘Er ġ,b,& 5 :с09ф) Ey 0,D,E5 “тг 0,0,6 5іп(ф) 


г%.б,Е с 





Ply Strain Strength Envelopes 








“ 


~ 
---.-.- е --е-- - -.- е.- -.-.-.....- 





Transverse Strain 





Longitudinal Strain 


—— 0 plv-Undamaged 

0 plv-matrix damaged 
----: 90 plv-Undamaged 
---- 90 plv-matrix damaged 

45 plv-Undamaged 

43 plv -matrix damaged 
-----. -45 plv-Undamaged 
---- -45 plv-matrix damaged 
—— Strain Relationship Line 
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Calculation of Degraded Laminate Stiffness: 
e first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 


stiffness calculated based on the ply damage history. 


plot:= fjale! 
ble! 
clel 
dle] 
counte 2 
niei 
while 1 


el ¡0 


51, –0 МРа 

Ex B.al,e, 
Ex Pß,bl+3.e, 
£ x DB,cl4-6.£; 
Ех В, 91+ 9,Е 5 
чап В) 


1 e 
El count 


£2 
count &l count 
] 
SSS See . ; М с М <. h N T а 
* statcount” KopNog+N nas tN p45 [No Q new], 1 N90 Qnewgj 5 ЊУ пао О пек, з EN pas Q newqi 4 ] 


B | 


ЕС —|| А 
Ylcount Stàfcount | | 


- m 


el count 


c2 
al count A erint count 


aleal+1 if еее xB.al,e, 


blebl+2 if ы, B,bl+3,e, 


clecl+1 if el ey B.clt+6.e 5 


die-di+1 if ER B.di+9,e, 


count count4- 

break if al --bl 4- cl 4-d1912 

break if al>3 

break if bI>3 

break if cl>3 

break if dl>3 
tempstore— augment a 

ER 


у] 
аџетеп! tempstore, = 
a 
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Ply Transverse Young's Moduli as a function of Strain 


Circumferential Ply Stiffness (Tensile) 


0 0 0 
0014 247801 16449 
0015 180.252 10.741 
0:016. 126253 707 
0.025 85233 2765 
0.171 94954 0554 


plot= 


Transverse Youngs Moduli (GPa) 











005 01 013 02 
Transverse Strain 


Transverse Stiffness Variations with Axial Location and Axial Crush 
istance: 


1. Setting up stations along the x' axis and axial crush distances at which E, is to be calculated 


п :=2000 ууу '=1..17 222:=1..п 


H 
ез еп, 
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2. Calculation of circumferential strain based on axial position and crush distance 


m prime; 9 temps H eff 1 


- nn — 


Е А ———— = 
0222 ууу R shell’ H etf ? tempus, 


3. Assignment of the appropriate Young's moduli to the correct axial position and crush distance based on the 
transverse Young's Moduli vs strain plot and circumferential strain. 


Eg :=|0-GPa if eg 2171 
ZZZ, yy) 222, ууу 

554GPa if 17l>e9 NEU 0 0 0 

S E 0.014 247801 16449 

2 GL - М-|0015 180252 10741 

| : SAN plot= = 

707GPa if Olé>e g,,, 2015 0.016 126.233 7.07 

> 014 0025 85.233 2.765 


10.741-GPa if 015>Е 922 ууу 0171 94954 0554 


16 449.GPa otherwise 


Flexural Strain and Stiffness Calculations, and Thickness Variations: 


1. Flexural Strain as a function of crush distance 0181 = 4.5729107 ? en 


2 +2 = 6 :=0,.001.. 2 АЈЕ. | 
ге 1 +20 5 , A hell t L 
с укб,ва):=2 PEF" 98. SAR һе ne 


^ Есен i 


Flexural Strain-O Orientation Plys 


0.08 


0.06 


0.04 


Flexural Strain 


0.02 





0 Бе | 
0 0 02 0.04 0.06 0.08 0.1 One 0.14 0.16 0.18 02 
Crush Distance 
---- РІУ1 
: Rives) 
Ss Plv 9 


..... Tensile Failure Strain for 0 Oriented Ply 


81 





Flexural Strain 


Flexural Strain 


01 


008 


006 


004 


002 


0.1 


008 


0.06 


004 


Flexural Strain-45/-45 Orientation Plys 


Lm me e o m mme 
—— td 


m 
"T 
аи 
e. 
lare 


vie ТУ 
ор ооо оо ти в ооо ооо оо ооо ооо оо фо фо фо фо оо оо ооо фо ооо ооо ооо оо фо оо фо ооо ооо ооо о ооо оо ооо ооо оо ооо о ооо ое ооо оо оо оо ооо 


, 
е 
е 


Am e <> <> <> <> <> <> <> <> <> <> <> <> => => <> <> => => => => => == 
‘ 


0.02 0.04 006 0.08 01 0.12 0.14 0.16 
Crush Distance 


Bin 

Ply 4 

Ply 6 

Ply 8 

Ply 10 

Tensile Failure Strain for 45/-45 Oriented Ply 
Compressive Failure Strain for 45/-45 Oriented Plv 


Flexural Strain-90 Orientation Plys 





0.18 





..... 


002 0.04 0.06 0.08 0.1 0.12 0.14 0.16 
Crush Distance 
Ply 3 
Ply 7 
Ply 11 


Tensile Failure Strain for 90 Oriented Ply 
Compressive Failure Strain for 90 Oriented Ply 


018 








2. Determination of flexural stiffness and thickness variations with crush distance. This a manual operation where 
the plys are failed progressively as crush distance is increased to 2H using the graphs above and the included ply 
strain limits. The thickness is determined by removing ply thicknesses as they fail. 








0 
0 
0 
0 
0 
0 
0 
0 
0 
: im Total plvs 
|, D РЕНЕ У E Q nexo [3 Re 2 | sso 
р x 
0 
р d star*=D star. | Е хр. = 9 star] | E E «£7 20 68439GPa 
0 
0 
0 
0 
0 
0 
0 
.002 164811 22h 
.003 16 4811 22h 
.005 17.3293 19-h 
.006 ` | 13.9983 15h 
.007 16.4377 13-h 
.008 18.3322 12h 
009 16 5229 11h 
ei 011 H eff ii 16.7947 К 10-h 
3 :=| 014 | —m E flexure *=| 17.0984 | GPa t=] Oh 
02 14 2239 8-h 
029 16 142 74 
041 17.2425 6h 
.068 18.3592 5h 
.093 12.799 4h 
H 18 4517 3-h 
™ 20 6843 2h 
2H 20 6843 2h 
m 
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Force vs Axial Crush Distance Determination: 


tL de >» Herp  EflexurenyR shell д Sv 


Eg Ex | [| ое 
„ ^рпт eff 
PPP, yyy “рр ^m эн 2 








rn x 
R shell Heft” ppp= 1 


Energy Absorption of the Structure: 


Епегоу( 611 ) *—linterp 6, P,611) 


УУУ 
Cumulatıve_Ener Ny = | ші Energy(x)d x Сштшайуе Епегеу оу 5) = 0 0219MJ 
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Appendix C 


Cylindrical Section Analysis 
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Ply Stiffness and Lay-up Inputs: 


1. Material System: E-Glass (1583)/Epoxy (7780) woven cloth 


2. Lay-up Configuration: [90,45,0,-45]3,90,45,45,90,[90,45,0,-45]3 


М 0:56 N 45156 

М 90 :=8 М p45 18 
3. Total Number of Plys Total_plys!= № о+М№ 245+ № р5+№ 90 
4. Ply Thickness h:=018.n 
5. Laminate Thickness t | = Тога _рјузћ 


6. Strain Applied to the Material (circumferential): 


B=atan(e,/e,) P :=89.5-deg 
Shear Strain: 150 


7. Number and Orientation of Ply Angles: 


Plys -=4 
0 
а.=1.. Руз 90 
Í 0:=| _ |-дев 

p 45 

c.=1..3-Plys - 
-45 
8. Engineering Constants: 
E 4 :73.MSI E , :7.55.MSI G m :=1.26-GPa 
= E m :=3 4-GPa 

E y :=2 7.MSI 

=> 

ум :=.5 

У Х'=.15 
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sequence .= 
a 4 


deg 








Component Parameters: 


1. Radius of Shell: 

2. Length of Shell: 

3. Wave Length of Elastic Buckling Wave: 
(shell only, no stiffener effect) 


4. Stiffener Parameters 
Undamaged Stiffness 


Number 
Area of Single Stiffener 
5. Open Section Parameters 


Reduction Parameter 


Ply Strength Parameter Inputs: 


1. Tensile Strength in Longitudinal Direction 


. Compressive Strength in Longitudinal Direction 


R shell :=3.5-ft 


L shell :=33-in 


H AR ѕһе L 


E carbon :=53.793-GPa 


N stiffener ‘8 


A stiffener ^ ( 5 3264n - ( 8 in— 2. 45-in))- 45-in 


:—62.5 deg 


Sí 


yes у =0.653 
л 





X :=39 1 KSI 


>< 


ргіте: 253 8-К51 


Y :=39.1-KSI 


. Compressive Strength in the Transverse Direction 


2 
3. Tensile Strength in Transverse Direction 
4 
5 


. Shear Strength 


Matrix and Fiber Damage Inputs 


E f star 01 : n y:7 316l п.-1 
Е := 15 :=316 F xy star 
паг‘ "> Ме = E 


Y prime 733 8KSI 


п 


:=9.41-KSI 


S 
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Damage Ply Calculations: 





Е у damaged F) -= 


. 1 1 
ку. = 4 
Е s damaged ) ~ Gf 


V 


m 
1+1 525 


m| — 
C 
а 


ито Ply Data: 
ndamaged ply, ply with matrix damage and ply with fiber damage. 


Е, Еу Es 
Е ух: Е E yy ‘=| © y_damaged © m_star E ss ‘=| Es damaged © m_ star 
EXE f star E y. damaged E f star E s damaged F f star 
x F xy star а E 
Ухх V xE m star FF xv. star :=| F ху яагЕ т 51аг n 
V XE f star F xy. star f star Š 554 





XX SX . 
primég prime’ Е Е 





ҮҮ =Y YY primeg *=Y prime т Y yyq 7 хха 


Rotation Matrices: 


ml, :—-cos 0 n, :—sin 
a a В = Ө; 
| e ә 5 
ES ns 5 . 
ml, a 2ml,n, 
J > 2 
=] п i. -2 n 
a a a 2ml, a 
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t3 
[29] 
to 


т! п. 2 ml, n, 4 ml na 
4 4 ^ 2 2 2 2 
n, т! 2 ml zn, 4 ml, an, 
2 2 | 2 2 1 4 4 1 2 2 
il. na m m m Tun mmi n 
Rotauonl, :— 5 
5 ee 24 2 M 2 ES m 31° 
ml, n4 ml n4 ami na [ mic un. | 
3 3 3 5 3 3 
ml, n -ml n4 ml пш ml n4 2[mi, n ml, ^] 
3 3 3 3 2 3 3 
| ml, na - ml na ml n Emi na 21 т! п mi, na || 
7 2 
ml, na 
Rotation? := 2 2 
а n, ml, 


Ply Stiffness Calculation: 


Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 


E хха у Ууд Е хха 
У 1-У ха У ууа 
4:5] У хха Ё ууд Е қ 
ШУ хха У УУа Iv хха ~ УУа 
0 0 E ss, 


9а 1,1 
One © тоха а | Qrotg a 2 O rota a E 
О ror, . :=Rotationl,.- = = 
rotq, a шо. Ms Q new, 477 Q rot à 5 O rota а > Qrotg a 5 
9% 33 So 5 O rota a 6 Ото а à 


Undamaged Elastic Stiffness of Laminate: 


NZ 1 
0 laminate ' 


a УНА Meet a =: | 
NotNo0tN nis tN pis [ 0 Q new; | +N 90 O new] 2 +N n5 © пем з +N p45 Q new, ] 


17575 5.256 -0.038 

a = = AF 

Q laminate=| 5256 17.726 -0.038 |°СРа SSO lanare E Eglass "3 
-0.038 -0.038 6.198 1,1 
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Effective Thickness Calculations: 


i = + A stiffener stiffener E carbon 
AE Eglass 


t eff = 0.035 t | =0.0128®m 
H eff! "A/R shell! eff H eff = 0 194m H=0.117%m 
IL, 
Number waves Tc Number waves = 2-159 
ені 


200 L shell 


Numb ec c H с e ЕАР 
umber waves eff Number waves H erp= 0.21m 


Flexural Stiffness: 
1. Indexing variables  ff:=1. Total plys4-! 
gg .=1.. Total_plvs 


2. Establishing an array of thickness positions for each ply 


SL 
PORT D 


4-- 
=. 





Zt 

3. Assigns the appropriate stiffness matrix based on the lay-up sequence 
One = Q new] | if sequence, 550 deg 

Q new] 2 if sequence, =90 deg 


Q new] ‚3 if sequence, 745 deg 


Q new| 4 otherwise 


Stress Space Strength Parameters: 
separate F must be caiculated tor each ply state 





1 ЕЕ ху_5їагу 5 
=. 
_ ^^ ^^ рите XX, XX prime, Y Ya YY prime, 
РЕЗЕ 
2 Xy. Star | 
ГЕ = зне n a У 
XX, XX prime, YYq YY prime; q ** primeg 
0 0 l 
55 2 

4 
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Strain Space Strength Parameters: 


1. Strain strength parameter for each ply state 


, 
E e F à 
al [Saa] рд бат 44 1.2 
О ат: + Е Го ото 
Be E - И с [ a 1,2] 
| | F В 
б F aq a af Oa 1 >] 2 9422 9% 2 9% nas 
и 
Gbaselq 7E P8 pa P a Pa an aa 
ME Q Q Q : ver ae Qn Wee 
811.21 Зал 94 2.2*[ 12] | л E 
"Faro 41,2 9422 “Ғаз2 22] 


Е я Е 
лы 


G = 
base}, Fe Q 


4 | а 1,27 15 Lu 


2. Calculating G matrices for all combinations of ply.state and angle: 


G base2| 
Rotation2, : 
G base2| 
G base?) 
Rotation2, · 
G basel ГІ 1 G base2> 
G basel 152 G base23 
Rotation]; : Rotation? . 
G baseli E G base23 
G basel} 3.3 С base2i 
G Rotation2, - G 
basel > 1.1 base2, 
G basel> 22 G base2, 
Rotation], - ' Rotation2, - 
Басе!» “| G pase2, 
- 1 „2 - 
G basel > 33 0 2 G base23 
otation2, - 
G basel; Ж G base23 
G 2trans = 
G basel; 33 | О ђазе2 


Rotationl, - E Rotation2, 


— 


—өә и 00 1 








G |trans' 


Rotation l : 


Когапоп |- · 


Rotationl., · 


Rotation], 


Rotation 1;: 


Rotation]; 


Rotation] 3! 


Rotationl 4 


Rotation] 4 


Rotation] 4 


С Базе! 
С вазе! ; 
С вазе! | 
С базе! | 
С Базе! | 
G basel} 
G basel> 
G basel> 


G basel> 


С Базе! > 
- 3 


G basel; 
G basel 3 
С вазе! 
G basel; 
G базе! | 
С Базе! | 


С базе 1 | 


G Базе! | 3 


G basel^ 
G basel^ 
G basel > 
G basel^ 
G basel 3 
G basel 3 
G basel 3 
G basel 3 
G Базе! | 
G basel 1 


G Базе! | 


G pasel 15 


С basel» 
G basel > 


С Базе! > 


С basel» 
=) 


С basel 3 


G basel 3 
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Rotation2, : 
С pase2, 


G base25 
Rotation2, - 
base2> 


G base2; 

Rotation2, - 
G pase? 

-3 


G Ба5е21 
Rotation2, - 4 
base2) 


G base2> 


Rotation2, - 
G базе2 


G base23 


Rotation2, - 
G base23 


1 


t3 








Rotation] 4 


basel; 12 


С Базе 15 


I 21,0 


G lirang | G I trans, 3 G l trans, 5 G 2rrans, | 
Gis G lirang , G | rans, 2 G lirang g G 2trans, 2 


G lirang 5 G гат 6 G Шап; | G 2trans. 3 


Formation of Strain Space Failure Envelopes: 


1. The following relationship is used to determine the failure envelope: G, уве) + бра 


2. After shifting to polar coordinates, the above equation can be solved for r using the quadratic equation: 


^ 
£ =F cosd Е.Г 5іпф Ar+Br— 1=0 


3. Calculation of coefficients for the quadratic formula 


A(0,b)'= б, узт) +2. б, р 9009) 9т0)+ 6 19959)" 


~ 
ar y 


Bela, т? G, sin( 6 )£ s4- G5 2 Ra б, | er Gb | ac 


x 
2,3 


~ 
С Ф.Б,= 5 = бь 3 a st Gp 3 3* s —1 


-B $.b,t, +B фЉе 2-4-А(ф.Ы)С ó. b.c, 
AS 


| Age au 
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ф :=0-дег, 1 -дег.. 360-дег 


Strain Space Failure Envelopes: 
cos В) 


у(х) :=хаж В) х х :20, е ..с09 8) 


Plv Strain Strength Envelopes 





^ 
` 
~ x 
D : 
eo Р = 
| ` — = 
` 
E ` 
~ a 
| 
` ^ 
% A 
| Se x 
| == FE ~ 
| An ZT Же ЗГА 
ES 6 "T CN ~ 
ra as O OO ааа” г .- af o.e -.-.-.--.- a corr a o n A A ow = 
БЕЗ ata aa nnn een een ; i 
` 
N % 
` : | 
7 eo ата” о => » 
У 7 с 
2 (eo 9 
` , | | 
` y | | | 
| 7 t+ ^ 
| ) t “ 
` ; | 
| ` 
` ; P | 
ql ; | 
/ 0 > 
` И | 
` | ; | 
A ; | 
, 
2 | 7 > ~ 
, 
-004 -0 03 -0.02.. 1-001 0 0.01 77 0.02 47 0.03 
^ y 
| ' a Nd 
“ t a : 
EI 7 
и А | 
о А 
УЖ / 
“ L] А 4 
~ X 2 
ns 
~ / 


, 
Е – - „u u -—— oom o mom om 
a e = e.c æa acam oe 





Transverse Strain 


Longitudinal Strain 


—— 0 plv-Undamaged 
: 0 plv-matrix damaged 
90 plv-Undamaged 
---- 90 ply-matrix damaged 
45 plv-Undamaged 
45 ply-matrix damaged 
T -45 ply-Undamaged 
---- -45 ply-matrix damaged 
—— Strain Relationship Line 
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Calculation of Degraded Laminate Stiffness: 
e first ply failure strain is first calculated. Subsequently, a progressive failure of plys is allowed and the 


stiffness calculated based on the ply damage history. 


plot:= fjale! 


ble] 
cle] 
dle! 
counte 2 
п] = | 
while l 
el ¡0 
61, =0-МРа 
су В.а1,е с 
En Bbl, Es 
el e min 
count ОВ. 
ex B.di+9e, 
ж? . 
& count =! count un B) 


A aM i qme 
stafcount NON GO+N nds tN pas 


-1 
-1 
Е певане || Á starcount de 


|» 0: © пема +% 90: Q newp] 2 +N pas Q newc] 3 +N pas: © пехар 4 | 


el count 
1е2 
Sl count] A statcount| “count 
aleal4-1 if el ounce x Balieg 


blebl+2 if el count? x B,bl+3,e, 


clecl+1 if 2! сое х В,С1+6,2 5 


dle-di+1 if clog х 8.91+9, 5 


counte count+ | 
break if al+bl4+cl+dl=12 
break if al>3 
break if bl>3 
break if cl>3 
break if dl>3 
tempstoree— augment лс 
МРа 
augment tempstore, 5 МІ 


‘GPa 
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Transverse Stiffness Variations with Axial Location and Axial Crush 
istance: 


1. Setting up stations along the x' axis and axial crush distances at which E, is to be calculated 


ууу.=1..18 п:=2000 zzz:=1..n 


002 
.003 
004 
.005 
006 
‚007 
| kc | 
mar, a 38 Нек 01 
Е „_| 013 | Нек 
о юра он 





.027 
.028 
.043 


115 
155 
234 


2. Calculation of circumferential strain based on axial position and crush distance 


EN prime, 9 temp, H eff 1 


£0 и 5 
ZZZ, yyy R shell H eff temp, 


3. Assignment of the appropriate Young's moduli to the correct axial position and crush distance based on the 
transverse Young's Moduli vs strain plot and circumferential strain. 
Ө ЈЕЈОФРа Кеб 2171 
ZZZ, yy) zzz,yyy 
469:GPa if .171>&g 2025 0 0 0 
222,У 


? 779.GPa if 025 = > 016 0014 245846 16.154 
2 (GPa 2 > 
б ж re 0.015 174.15 10.133 


0 
6 094.GPa if Olé>e 0.2 ууу 015 0.016 113.553 6.094 
0.025 90.216 2779 
0.171 80.302 0.469 


IV 
= 
> 


10133GPa if OlS>eg а 


16.154-GPa otherwise 
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Flexural Strain and Stiffness Calculations, and Thickness Variations: 





1. Flexural Strain as a function of crush distance 


2 +2 52 5 :=0,.001. 24/Е 1 
ee+ 17 “gg б e >> shell * L 
= A19 A[R слеп + 


Е xfí 9,88) Fre 


Flexural Strain-O Orientation Plys 





= 
= 
© 
E 
? 
ES 
1232424..525221222.7 
0 0.05 0.1 0.15 02 
Crush Distance 
— Ph | 
use Ply 5 
= РІУ 9 
..... Tensile Failure Strain for 0 Oriented Ply 
<----- Compressive Failure Strain for O Oriented Ply 
, Flexural Strain-90 Orientation Plys 
0.2 
0.15 
E 
© 
© 
= Öl 
= 
ж 
ED 
0.05 
0 





Crush Distance 


---- Ply 3 

Ply 7 
T Ply 11 
42... Tensile Failure Strain for 90 Oriented Ply 
-=- Compressive Failure Strain for 90 Oriented Ply 
сзсз Tenstle Failure Strain for 45/-45 Oriented Ply 
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Flexural Stratn-45/-45 Orientation Plvs 


Flexural Strain 


0 005 0.1 0.15 
Crush Distance 


— Ply2 
Ply 4 
---=-- РІУ 6 
Sc РІУ 8 
== Рју 10 
4... Tensile Failure Strain for 45/-45 Oriented Ply 
= = - Compressive Failure Strain for 45/-45 Oriented Ply 
сәсе (ШӘ8Ы 
—— trace 9 


2. Determination of flexural stiffness and thickness variations 
with crush distance. This a manual operation where the plys are 
failed progressively as crush distance is increased to 2H using 
the graphs above and the included ply strain limits. The 
thickness is determined by removing ply thicknesses as they fail. 


Total plvs 
12 


3 3 
3 2 TS [E xxl y] шт; Б 
(24) хх = | 


Lal — 


D star := 





= -1 
d Star D star 


Курт Ч заг | 


E „f= 11.38740GPa 
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16.4023 
16.4023 
16.7798 
15.9182 
17.3134 
16.0961 
16.4632 
16.0174 
16.5986 
19 1109 
16 705 

17 9956 
17.8391 
13.7589 
16 6569 
17.9459 
231851 
11.3874 


-GPa 


55: 


е, С) fe) te) С с: сз со со соо Лы «Ооо о Оо оо о ос осо 








Force vs Axial Crush Distance Determination: 


10 


3 
tet Hef F Пехигеууу К shell’ ! f, 


Ее Bg 
К слеп Не ppp= 1 


Energy Absorption of the Structure: 


Епегеу( 611) :-ішшегр(6,Р,611) 


EX 1 p 
ррр.ууу ^ Püm€ppp 7m x 3 
3. H eff 


€ Е, . 


| буу, 
= Energy(x)dx Ситшапуе Епегеу = 0 244°М] 
om rows à ) 


у oec У 
‚Cumulatı e Energy y: Cumulati e Energy. _ 


ж ws Е о": 
instant_Energy, уу if vv], Cumulatv e Energy... S 
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Appendix D 


Stiffener Analysis 
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Ply Stiffness and Lay-up Inputs: 


1. Material System: Carbon, woven cloth (6K 5HS satin weave IM7) 


2. Lay-up Configuration: [90,45,0,-45]4[90,45,0,-45], 


N 0 :=8 М п45 :=8 

N 99:78 М5: 
3. Total Number of Plys Total_plys'= N 9+N p4s+N p35+N 90 
4. Ply Thickness ћ:= 0135 
5. Laminate Thickness t | :=Total_plysh 


6. Strain Applied to the Material (circumferential): 


B=atan(e,/e,) 8 :=89.5-дев 
Shear Strain: 6150 
7. Number and Orientation of Ply Angles: 


Рју5 :=4 


0 
a :=1.. Plys ._| 90 
0:-|  |deg 
:=1..3.Ріуѕ ps 
c:=1..3-Ph 245 
8. Engineering Constants: 
E , ‘=11.6-MSI E,:737MSI G m :=1.26-GPa 
es E m :=3 4GPa 
E y :=11.6-MSI 
MAS 
ут .=.5 
v x. 7.124 
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Component Parameters: 


1. Radius of Shell: 
2. Length of Shell: 


3. Wave Length of Elastic Buckling Wave: 


4. Number of Buckling Waves Along Shell 


Ply Strength Parameter Inputs: 


1. Tensile Strength in Longitudinal Direction 


. Compressive Strength in Longitudinal Direction 


К пей 2350 


L shell 1544710 


H AIR shell! L 


¿_L shell 
Number waves ж 


X :=39 1-KSI 
X prime =33 8KSI 


У :=39 1 К5! 


. Compressive Strength in the Transverse Direction 


2 
3. Tensile Strength in Transverse Direction 
4 
5 


. Shear Strength 


Y prime 733 &:KSI 


Matrix and Fiber Damage Inputs 


Е є star = 01 пу:= 5161 n:=1 


Е то чаг' = 15 n s:=316 


Damage Ply Calculations: 








E y_damaged F).= 


E s damaged F)':7 


F xy star 


$ .=9 41.К $1 
= 
“т “т 
p Al rer 
ERE = 
p ES ЕТ 
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PIRE Ply Data: 
ndamaged ply, ply with matrix damage and ply with fiber damage. 


Ex 
Ехх:=| Fx 


EXE f star 


Vy 


У хх У ХЕ т star 
V x E f star 


q md xm primeg SP prime 


Rotation Matrices: 


ml, +=c0s 9, па яп 9. 
D 
ml, na 2ml.n, 
uu = mi.” -2m 
а а m тати, 
l | в 2 
min, min, т zen 
4 4 
ml n 
a a 
n + ml : 
a a 
5 = a 
ml, ^. n. ті .n 
a 
Rotation], ‘= 
а 222 2 
ml n ml -n 
a 
3 
ml, т, а n, 
5 3 
ml -n - ml : 
a a 
| 2 2 
т n, 
ion? = 2 2 
Rotation2, ; п mit 


кә 


~ 


у 


Еуу:=| Еу damaged Fm star 


E y_damaged E f_ star 


F xy. star 


FF xy. star? 7| F xy_starE m_star 


F xy. starE f star 


Es 
E ss :=| E s_damaged E m_star 


E S damaged E f star 





7 2 q 
XX rim =х ri e 
primeg ргіте -р- 


e 
“2 
о 





Ее 


о 
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Ply Stiffness Calculation: 


Calculates the zero direction stiffness matrix for the various ply states (undamaged, matrix damaged, and 
fiber damaged) 


Е VIC ER 
х заха o 
1-- У ха ~ Ууа |-у XX AT YYq 
к= ОУ ү Ел. Е 
7, “Буға Ууд 4 
|- > У Ууд Г-У хуу yq 
0 0 E 
554 


ати 
Oe Orotga y Qreqa, Qro, 
Q rot, , /"Rotationl,- dli pal 
TOig a Ё p (5 Q new, қ. rots 3 Оша ә а 6 
9% E Q rotg a Р Что а 2 Что, а E 


Undamaged Elastic Stiffness of Laminate: 


Q laminate =N EN FN ms FN pas |, 0: Qnew, | +N до: Q new) > +N 145 О new) 5 + № р45: Q new, , ] 


= 





64 715 26 586 63220715 
жү 27 -i ү ым 
О laminate 7 | 26 586 64 715 –6.322%107 !> јоРа  5:=0 јатпаке S ES 
632210715 -632240715 19064 
| =53.793°GPa 
35.2 
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Appendix E 


Forward Hemi-head Pressure Hull 


Analysis 
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Analysis of HY-80 Forward Pressure 
Hull Hemihead 


Stress Strain Relationship for HY-80 


0 0 

‚002567 77 

.0027 78 

E- 0029 ик 79 KSI 

0031 79.5 

.0046 80 

.01 80.5 

L .02 | 81 
rows E) 


х= (0, curve(x) :=linterp(e,o,x) 


7100 ""rwss) 


Circumferential Ply Stiffness (Tensile) 


600 














Stress (MPa) 








0 0.002 0.004 0 006 0.008 001 0.012 0.014 0.016 0.018 0.02 
Transverse Strain 
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Material Flow Stress Determination: 


1 E rows =) 


Qc curve(x)dx 





€ 
TOWSÍ £ ) 0 


Structural Properties: 


Е 7 
tL тен | 


bbb :=1..rows(ó) 


Crushing Force Prediction: 


tal us 
tal — 


Pob = т 42-5 с б, 


Energy Absorption: 
Energy(611) ;z linterp( 6, P,611) 


“бы, 


Cumulative Energy, :7 


Crushing Force for one Wave Length 


10 
E 
"= 
> 
2 5 
< 
en 
2 
0 
0 0.1 02 0.3 


с c 7.517.499 «MPa 


Energy(x)dx 


0.4 


Axial Crush Distance (m) 
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0.5 


0.6 


0.028 
0.1 
0.241 
0242 
0.245 
0.247 
0.248 
025 
0252 
0.254 


0.258 | 


0.264 
0.268 
(277 
0.289 
p.29] 
0.318 
0.349 
0.447 
0.519 
0.66 








Cumulative_Energy = 9MJ 
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